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Summary 

The merits and limitations of the optimality criteria (OC) 
method for the minimum weight design of structures subjected 
to multiple load conditions under stress, displacement, and 
frequency constraints were investigated by examining several 
numerical examples. The examples were solved utilizing the 
OC Design Code that was developed for this purpose at NASA 
Lewis Research Center. This OC code incorporates OC 
methods available in the literature with generalizations for 
stress constraints, fully utilized design concepts, and hybrid 
methods that combine both techniques. It includes multiple 
choices for Lagrange multiplier and design variable update 
methods, design strategies for several constraint types, variable 
linking, displacement and integrated force method analyzers, 
and analytical and numerical sensitivities. On the basis of the 
examples solved, this method was found to be satisfactory for 
problems with few active constraints or with small numbers of 
design variables. The derivable OC method without stress 
constraints was found satisfactory even for large structural 
systems. For problems with large numbers of behavior 
constraints and design variables, the method appears to follow 
a subset of active constraints that can result in a heavier 
design. The computational efficiency of OC methods appears 
to be similar to some mathematical programming techniques. 

I. Introduction 

Three procedures are currently available for automated 
structural design: (1) the fully utilized design (FUD) (refs. 1 
and 2), (2) the optimality criteria method (OC) (refs. 3 to 5), 
and (3) the mathematical programming techniques of opera- 
tions research (MP) (refs. 6 to 9). The FUD, because of its 
simplicity, is popular in industry despite the availability of the 
other two design methodologies. FUD may, however, exhibit 
some limitations, especially for stiffness and dynamic 
constraints (ref. 6). To alleviate deficiencies associated with 
FUD, the structural design problem was formulated and solved 
during the 1960’s as a nonlinear mathematical programming 
optimization problem. Structural mass was used as a typical 
objective and failure modes as the constraints (ref. 6). This last 
approach is conceptually sound and mathematically elegant, 


but its solution can be computationally prohibitive. At about 
the same time, a relatively simpler design technique, termed 
the optimality criteria method and based on a Lagrange 
multiplier approach, was introduced to solve certain types of 
structural optimization problems (refs. 3 to 5). 

Theoretical aspects of the optimality criteria method have 
already been given (refs. 3 and 4). In this report the merits 
and limitations of this method, which may provide a practical 
design tool, are considered by examining several numerical 
examples. The examples are solved using the optimality cri- 
teria design code that was developed to assess the perfor- 
mance of this design technique. 

The optimization section of the OC code contains numer- 
ous choices for optimality criteria update formulas, including 
most of the formulas in references 3 and 4, along with fully 
utilized design rules. The optimality criteria method for 
structural optimization was originally derived (refs. 3 and 4) 
for displacement buckling and frequency constraints. Stress 
and frequency constraints, however, are also included in the 
OC code, because such constraints are relevant to structural 
design. The analysis segment of the code includes a choice of 
three finite element analysis methods: (1) the displacement (or 
stiffness) method, (2) the integrated force method, and (3) a 
simplification of the integrated force method. The OC code 
has options to calculate design sensitivities of the behavior 
constraints either in closed form or using numerical differ- 
entiations. To compare results obtained by the optimality 
criteria methods with those of other structural optimizers, we 
integrated the OC code into the Comparative Evaluation Test 
Bed of Optimization and Analysis Routines for the Design of 
Structures (CometBoards), an optimal structural design test 
bed under development. CometBoards offers several choices 
for optimizers and analyzers. 

The subject matter of this report is presented in the sub- 
sequent five sections. In section II, the theoretical aspects of 
the optimality criteria method are introduced. In section El, the 
structure of the OC code software is briefly discussed. In 
section IV, several numerical examples are provided. The 
issues associated with the optimality criteria method, as well as 
merits and limitations of the three analysis methods, are 
examined and illustrated in section V. The conclusions are 
given in section VI. Lists of symbols and of acronyms and 
initialisms are provided in appendixes A and B, respectively. 
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II. Theoretical Bases of Design Methods 

The three automated design procedures — (1) optimization 
using mathematical programming techniques, (2) the 
optimality criteria method, and (3) fully utilized design 
concepts — are briefly examined and compared. 

A. Mathematical Programming Techniques 

In mathematical programming techniques, structural 
design is cast as the following optimization problem: 

Find the n design variables % within prescribed upper 
and lower bounds [xt^ -Xi ~xF S > 1 = that 

make the scalar objective function f ( % ) (here, weight) an 
extremum (minimum) subject to a set of m inequality 
constraints, usually defined as the failure modes of the 
design problem: 

gj(X)~ 0 (7 = 1,2,..., m) ( 1 ) 


in, say, K design updates. At each iteration the design is 
updated by calculating two quantities: a direction ipand a 
step length a. The optimal design process, utilizing the direc- 
tion and associated step length, can be symbolized as 

Jt=i 

where (p k is the direction vector at the £th iteration and a \ is 
the step length along the direction vector. 

The direction vector at the £th iteration is generated from 
the gradients of the objective function and the active con- 
straint subset following one of the many available direction- 
generation algorithms (refs. 10 and 11). Along the direction 
vector Jp k , a one-dimensional search is carried out to obtain 
the step length a again utilizing one of several available 
procedures (refs. 12 and 13). The updated design is then 
checked against one or more stop criteria (ref. 1) until it con- 
verges. The details of the nonlinear mathematical program- 
ming techniques, well documented in the literature, are not 
elaborated here (refs. 10 to 15). 


The constraints for structural design applications are 
typically nonlinear in the variables % , thus it becomes a non- 
linear programming problem. Note, here, that equality con- 
straints could also be included. 

This report considers stress, displacement, and frequency 
constraints (gj) under multiple load conditions. For each load 
condition, the stress constraints are specified by 


B* Design Cast in a Lagrange Multiplier Formulation 

The optimality criteria method is an alternative design tool 
for solving the optimization problem given by equations (1) 
to (4). This method can be considered to be a variant of the 
Lagrange multiplier approach, which has been specialized for 
structural design applications. We consider next the Lagrange 
multiplier method for solving the optimization problem. 

The Lagrange multiplier approach adjoins the constraints 
to the objective function, and the Lagrangian is formed: 


where Oj is the design stress for the jth element and Ojq is 
the permissible stress for the jth element. For each load 
condition, the displacement constraints are specified by 


Sjs+j 



1 . 0*0 


(3) 


where uj is the jth displacement component, ujq is the dis- 
placement limitation for the jth displacement component, and 
js is the total number of stress constraints. Constraints on 
frequencies are specified by 


g f = 


fnO _) 2 

. fn j 


- 1 . 0<0 


(4) 


2(?.*)-/Cf)+ 2 <6) 

active set 

where the superscript asterisk indicates those constraints 
within the active set. 

The passive constraints do not influence the design, and 
associated Lagrange multipliers are zero. The independent 
variables associated with the Lagrangian function (eq. (6)) are 
the design variables and multipliers associated with the active 
constraints. These two sets of unknowns, the design % and the 
multipliers X , in principle can be obtained from the stationary 
conditions of the Lagrangian, =$?(£, A ). The stationary 
condition of the Lagrangian with respect to the design 
variables yields 

wh 2^ v »; ®-° <7> 

active set 


where f n represents natural frequencies of the structure and 
f n0 the limitations on these frequencies. 

The optimal design £ opt in a mathematical programming 
technique is obtained iteratively from an initial design 


Likewise, the stationary condition with respect to the mul- 
tipliers yields the active constraints 

g * [x] = 0 (g/ within the active set) (8) 
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Simultaneous solution of equations (7) and (8) would 
yield the solution to the optimization problem. For structural 
design problems, this direct approach can be computationally 
intensive and may suffer from numerical instability. 

The optimality criteria method bypasses direct Lagrangian 
solution. Instead, the intrinsic nature of the design problem 
has been exploited, and the process has produced several pro- 
cedures. The derivation of an optimality criteria method for 
minimum weight optimization of structures in the context of 
displacement constraints is illustrated in appendix C. 

Before the optimality criteria methods are presented, the 
relation between the number of design variables and the 
number of active constraints, which plays an important role 
in the optimization process, is examined. The three possible 
relationships between the number of design variables (NDV) 
and the number of active constraints (NAC) of the optimal 
solution are illustrated by considering a three-bar truss 
(fig. 1). The bar areas are taken to be the three design vari- 
ables. The truss is subjected to two load conditions, which 
give rise to 6 stress, 4 displacement, and 1 frequency con- 
straint, or a total of 11 inequality constraints. 

C. Relationship Between Number of Active Constraints 
and Number of Design Variables 

The three cases, (1) NAC greater than NDV, (2) NAC equal 
to NDV, and (3) NAC less than NDV, are examined separately. 
In the discussion of these relationships that follows, functional 
independence among the active constraints is assumed. Note, 
however, that functional dependence among the behavior 
constraints of structural systems can occur and may require 
special consideration as elaborated in reference 16. 

L When the number of active constraints is greater than 
the number of design variables (NAC > NDV ). — Consider the 
design space of a problem with n design variables and n + v 
active constraints at the optimum. The optimal solution, a 
point in that design space, can be specified as the intersection 
of any n of the ( n~ l)-dimensiona! surfaces defined by the 


active constraints, assuming that the active constraints are 
locally well behaved and functionally independent. In a more 
general sense, if there exists a subset, of any size, of active 
constraints the intersection of whose surfaces form the optimal 
solution point, then there exists a set of, at most, n of those 
surfaces whose intersection is sufficient to define the optimal 
point. The remaining v constraints, which pass through the 
optimal point, may be termed “follower constraints ” Such an 
occurrence is illustrated in Figure 2(a). 

For a three-bar truss with three design variables, consider 
an optimal design that has five active, functionally indepen- 
dent constraints, g 3, £5, £5, £9, and gn, as depicted in 
figure 2(a). Any set of three active constraints is sufficient to 
establish the optimal design (for example £3, £5, and £$, as in 
fig. 2(b)). The two remaining active constraints (here, £9 and 
£n) are the follower constraints and need not be considered in 
optimization calculations. In summary, from geometrical con- 
siderations, the inclusion of a maximum of n active constraints 
in the optimization process is sufficient to generate the optimal 
design vector of dimension n, By considering only as many 
active constraints as there are design variables, we reduce the 
complexity of a given optimal design problem. 

2. When the number of active constraints equals the 
number of design variables ( NAC = NDV ), — The solution of 
an optimal design problem in which there are n design 
variables and the same number of active constraints (i.e., 
NDV = NAC = n) is, by definition, a fully utilized design. 
The stationary conditions of the Lagrangian 

^{xX)= f{x)+Y^x jgj {x) 

j = 1 

for this situation yield the following equations: 
n 

V /(*) + X V V &*j(x) = 0 (9a) 



Sj(x) = ° (j = l,2,...,n) (9b) 

Equations (9a) and (9b) together represent 2 n equations in 
In unknowns: n design variables (X/l and n Lagrange multi- 
pliers {A,/}. The stationary condition of the Lagrangian with 
respect to the des gn variables (eq. (9a)) yields n equations in 
2 n unknowns (i.<\, n Lagrange multipliers and n design 
variables). The stationary condition of the Lagrangian with re- 
spect to the multipliers reproduces the n constraint equations in 
n unknowns (eq. (9b)), which are independent of the Lagrange 
multipliers. The set of n constraint equations (eq.(9b)) alone is 
sufficient to generate the n design variables, if we assume that 
they are functionally independent. Once the design variables 
are known, equation (9a) can be used to calculate the Lagrange 
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truss. 



(b) Three active constraints, which are sufficient to determine the 
optimal point, and two follower constraints. 

Figure 2. — Optimum design point. 


multipliers. It is important to note that the objective function 
does not participate in this calculation of the design. In the 
case of the three-bar truss example, if exactly three constraints 
are active, then simultaneously solving the three constraint 
equations (eq. (9a)) would produce both the optimal and the 
fully utilized design. Equation (9a) can then be used to calcu- 
late Lagrange multipliers, if desired. 

In summary, when an optimal design has functionally inde- 
pendent active constraints that equal or exceed the number of 
design variables, then the design represents both a fully uti- 
lized design and an optimal design. 

3 . When the number of active constraints is less than the 
number of design variables (NAC< NDV ). — If the number of 
active constraints in the optimal design is less than the number 
of design variables, the design is not fully utilized. Some fea- 
tures of this case can be illustrated by the three-bar truss. Let 
us assume that, at optimum, two constraints, g\ and g 2 , are 
active. The stationary conditions of the Lagrangian yield the 
following equations: 

w{x)+^s*{x)+^yg2(x)= 0 ( 10a ) 

S/(*) = ° (j = 1- 2, •••.«) (10b) 

Equations (10a) and (10b) represent five equations in five 
unknowns (three design variables, X2> ancl and two 
Lagrange multipliers, X\ and X 2 ). The stationary condition of 
the Lagrangian with respect to the design variables (eq. (10a)) 
represents three equations in five unknowns (the two Lagrange 
multipliers and the three design variables). The two constraint 
equations given in equation (10b) represent two equations in 
three unknowns (the three design variables). Although these 
last two equations are independent of the Lagrange multipliers, 
they are insufficient in quantity to solve for the three design 
variables. Thus, all five equations (eqs. (10a) and (10b) are 
coupled in the design variables and the multipliers, and the 
entire set must be solved simultaneously to generate the solu- 
tion. The gradient of the objective function and the Lagrange 
multipliers do participate in the calculation of the design vari- 
ables because of the intrinsic coupling between equations 
(10a) and (10b). In other words, only when the number of 
active constraints is fewer than the number of design variables 
do both the constraints and the objective function participate in 
the optimization, this situation is the one most frequently 
encountered. 

D. Optimality Criteria Method 

The optimality criteria method provides several procedures to 
iteratively update both the design variables and the Lagrange 
multipliers. The update rules described below include modifica- 
tions and generalizations of the original formulas given in 
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references 3 and 4 for stress, displacement, and frequency con- 
straints. The derivation of one rule is illustrated in appendix C. 

The active constraint set is estimated through a constraint 
thickness parameter that defines a finite interval [— fy, 0] 
within which all constraints are considered active. This con- 
straint thickness becomes progressively tighter at each suc- 
cessive iteration, ultimately reducing to zero. The constraint 
thickness update rule used is 

t k = r k ~ x t 0 (k = 1, 2 , K) ( 11 ) 


i*+i - l k 
A ja “ A ja 


f r \ a Po 

Z E 

r* 

k l j°j 


( 13 ) 


The optimality criteria formulas (eqs. (12) and (13)), do not 
use weight in calculating multipliers. 

c. The Lagrange inverse forms: The multipliers in these 
methods are obtained from equation (7) by premultiplying it 
with the transpose of the gradient matrix Vg and then 
solving for the Lagrange multipliers: 


where f* is the constraint thickness at the kih iteration, to is 
the initial prescribed constraint thickness, t(< 1) is the factor 
by which the constraint thickness is reduced at each iteration, 
and K represents the number of iterations required to reach a 
solution. All violated constraints whose values are greater 
than zero are also included in the “active” constraint set. 

L Lagrange multiplier update formulas . — Utilizing 
equations (eq. (7) and (8)), references 3 and 4 give several 
update rules for the Lagrange multipliers, two of which are 
included here. Three additional update methods are also 
discussed. Refer to appendix D for a summary of the update 
formulas. 

a. The linear form: The linear form of the update formula 
for the Lagrange multipliers is obtained by manipulating the 
constraint equation and has the following form: 

X k ; a x =X%[l.O + cc k p 0 (C ja -C* a )] (12) 

where X k a is the value of the Lagrange multiplier associated 
with thcjath “active” constraint at the kih iteration, Cj Q is the 
actual value of the displacement at a particular node, or the 
stress in a particular element, and Cj a is the corresponding 
permissible value of that displacement or stress (i.e., 
Cj a - (gj a + l)C} fl ). For frequency constraints, Cj a is taken to 
be the square of the inverse of the natural frequency, 1 //„ 2 , 
and C*j a is l/f n o 2 - Also in equation (12), p$ is the initial value 
of an update parameter, and a is an acceleration parameter 
used to modify that update parameter. The value of po is typi- 
cally 0.5, and a is often taken as 1 .0. 

Only those Lagrange multipliers are updated which are 
associated with active constraints. Active constraints are 
defined on the basis of the thickness parameter given by 
equation (11) (i.e., Xj a , where ja is such that gj Q ~0c)- 
However, all Lagrange multipliers may, in some cases, be 
updated, regardless of the corresponding constraint value. 
Initial values for the Lagrange multipliers are required to begin 
iterations. 

b . The exponential form: The exponential form is also 
obtained from the constraint equation and has the following 
form: 


r^lV]] 1 [vrfv/ (i4) 


where A* is the vector of Lagrange ljiultipliers associated 
with the “active” constraints, V/ is the gradient of the 
objective function, and the superscript asterisk represents the 
fact that only active constraints are included in the sensitivity 
matrix. If the number of active constraints m exceeds the 
number of design variables n, then the number of columns m 
in Vg*will be larger than the number of rows (n; m > n), 
making the m X m coefficient matrix [[ Vg*] r [ Vg * ]] 
singular. As mentioned earlier, however, a maximum of n 
active constraints is sufficient to establish an optimal solution 
in an n-dimensional design space. Furthermore, except for 
functional dependence among active constraints, restricting 
attention to any n of the active constraints at optimum should 
yield the optimal solution. With this Lagrange inverse form 
of the multiplier update method, no more than n active 
constraints are considered at any given time. 

Two more Lagrange multiplier update formulas can be 
derived — the unrestricted and the diagonalized inverse forms. 
In the unrestricted Lagrange inverse form, equation (14) is 
used, but the number of active constraints considered is 
allowed to exceed the number of design variables. This form 
is included to examine the difficulties that may arise from 
singularities whenever the number of active constraints 
exceeds the number of design variables. 

The diagonalized Lagrange inverse form is obtained by 
simplifying equation (14) (i.e., only diagonal entries in the 
matrix [Vg r [Vg *] are taken). The coupling terms in the 
sensitivity matrix are neglected. This update rule becomes 


Xj a — 


Vgja 


(15) 


2. Design variable update formulas . — The optimality 
criteria update formulas utilize a vector 25 of dimension n. 
The vector D is calculated from the stationary condition 
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(eq. (9)) by scaling, component by component, with respect 
to the gradient of the objective function and has the following 
form: 

^i a 8 J a )/ 

Dj = ~ ja — (i =1,2 n) (16) 


where j a represents the ith component of the 

summation of the product of the gradients of the active 
constraints and the Lagrange multipliers, and V ft is the ith 
component of the gradient of the objective function. The 
stationary condition of the Lagrangian with respect to the 
design variables (eq. (9)), in terms of the vector 15, can be 
represented as 


D f = LO (17) 


Design updates that are based on equation (17) (refs. 3 
and 4) are given next. Refer to appendix D for a summary of 
the design update formulas. 

a. The exponential form: The exponential form of the 
design variable update method is as follows: 


1 _ v k n \/(p k q 0 ) 
Jii Ai L'i 


(18) 


where 40 and (i are acceleration parameters. Typical values of 
#0 and fi are 2.0 and 1 . 0 , respectively. 

b . The linearized form: The linearized form of the design 
variable update method can be considered to be a truncated 
Taylor expansion of the exponential form (eq. (18)) and is 
given by 


v *+l _ 

41 Ai 


1 ° + (A- 1 0 ) 

P % 


(19) 


c. The reciprocal form : A Taylor expansion of the inverse 
of the exponential form (eq. (18)) yields the reciprocal form: 


Xi 


*+i 


Xi 


1 . 0 -- 




(A -1.0) 


( 20 ) 


d. The melange form: These rather complicated formulas 
to update the design variables are an attempt to combine 
influences of several factors simultaneously: active constraint 
types, the rescaling vector Z5, and factors to stabilize con- 
vergence. The design variable is updated in two steps. The 
initial step, represented by k + / 2 , is 



( 21 a) 


The final step is 


v *+i _ 

Al 


MX 


*+ 1/2 


+ Xi 0-°+g ( .s)] 

yk+ 1/2 


if Xie{xi\gis*-tk} 

Xi « {x\8is^~ t k} 
(21b) 


In equation (21a), the coefficients, a , b , c, and d depend on 
active constraint types (stress, displacement, or frequency), as 
well as certain prescribed parameters. The coefficients, how- j 
ever, always sum to unity. Here, g ^ and ^ fl / are constraint 
values for specific “active” displacement and frequency con- - 
straints, respectively; njad and njaf are the number of each of 1 
these types of constraints in the “active” set; and +U2 is an 
intermediate design value that may be modified within a given 
iteration. Only active constraints are considered in this update | 
formula. 

The second term in equation (21a) shows a sum over the 
displacement constraints included in the “active” set. If there 
are none, this term is dropped. The factor, 1 + g, is the scaling 
factor that will move the design just inside the feasible 
domain with respect to any one displacement constraint. Note 
that the largest such factor is the one used in the fully utilized 
design methodology described later. Although stress con- 
straints are closely linked to the cross-sectional area of the 
element corresponding to each constraint, displacements are 
global responses and can be significantly influenced by any 
or all design variables. The third term shows a sum over all 
frequency constraints included in the “active” set and also 
represents a global response. The last term uses the rescaling 
vector 15, with an exponent of one, and the first term retains 
the old (prior iterate) value of the design variable. 

Any design variable associated with an element whose 
corresponding stress constraint is included in the “active” set 
is further modified by taking the arithmetic mean of the value 
obtained in equation ( 21 a) with the value of that design vari- 
able that would be obtained by a fully stressed (fully utilized) 
design technique (see eq. (21b)). Here, gj s is the value of the 
stress constraint corresponding to the ith element. Note that 
only one element is influenced by any one active stress con- 
straint. All design variables not associated with any such ele- 
ment retain their values given in equation ( 21 a). 

3. Hybrid design variable update methods . — The hybrid 
methods represent an augmentation of the optimality criteria 
design update rules to the fully utilized design concept. These 
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three methods were devised in an attempt to retain the ben- 
efits of the fully stressed design for stresses with the 
optimality criteria method for displacements and frequencies. 
For these methods, the rescaling vector D , given by equa- 
tion (16), is calculated for active displacement and active 
frequency constraints only. To initiate the design variable 
update, each method uses a corresponding optimality criteria 
update formula (eq. (18), (19), or (20)) to obtain an interme- 
diate design. Next, an alternative intermediate design is 
obtained with a fully stressed design for the stress 
constraints. 

Xi +l/2 =xH l0 + 8is) ( 22 ) 

Whenever more than one stress constraint is associated with 
any given design variable (e.g., with multiple load condi- 
tions, or linking of design variables), the largest of the scaling 
factors (1 + g) is used. Finally, the two intermediate designs 
are compared, one design variable at a time, and the largest is 
chosen. These three design variable update methods are 
termed (1) the hybrid exponential form, which incorporates 
the exponential form of the optimality criteria design update 
formula (eq. (18)); (2) the hybrid linearized form, which uses 
the OC linearized form (eq. (19)); and (3) the hybrid recipro- 
cal form, which uses the OC reciprocal form (eq. (20)). 

E. Fully Utilized Design Approach 

The fully utilized design is obtained in two steps. First, a 
fully stressed design, designated ;£ a ’ opt , is obtained by using 
the stress ratio technique. The design is then rescaled for vio- 
lated displacement and/or frequency constraints, if any. The 
two steps are — 

L Fully stressed design . — For a fully stressed design, the 
design variables are updated for stress constraints only, as 

Zi MX =Xi**<* (*• = 1.2. -.») (23) 

where xf' k the ith component of the fully stressed design 
at the A;th iteration, and the scaling factor for the ith 
design variable is 

R ma x(g Ut g a °u) (24) 

where G\ /, <72/, <7u represent the stress values for each 

element associated with Xi under each load condition. Note 
that linking of design variables can cause more than one ele- 
ment to be associated with a given design variable for a given 
load condition (see section III.D). The fully stressed design 
^cr.opt j s known converge in a few cycles, irrespective of 
the number of design variables (ref. 1). 


2. Fully utilized design . — The fully utilized design for 
simultaneous stress, displacement, and frequency constraints 
is obtained in one additional step by scaling the fully stressed 
design to satisfy the maximum violated displacement and fre- 
quency constraint, if any: 

r^r-oP'O + Smax) (25) 

where g max is the value of the largest violated displacement 
or frequency constraint. 

III. Optimality Criteria Computer Code 

The OC computer code is composed of three modules: 
(1) the optimization module, (2) the analysis module, and 
(3) the interface module. The code was developed on the 
basis of existing interfaces to analysis routines and was incor- 
porated into CometBoards, an optimal structural design test 
bed. 

A. Interface Module 

The interface module reads the problem specification, 
such as the geometry, material properties, design limitations, 
and optimization parameters from free-format keyword- 
based data files. 

Interactive input of some of the optimization parameters is 
an additional optional feature. This module also initializes 
variables, calls the optimizer, and prints out the final results. 

B. Optimization Module 

The optimization module of the OC code was designed to 
provide considerable flexibility, such as multiple choices of 
both Lagrange multiplier and design variable update methods 
and parameters, as well as various strategies for approaching 
problems with several constraint types. 

The core iteration loop of this module proceeds as 
follows. The thickness for each constraint type (stress, dis- 
placement, and frequency) is modified by a constant factor, 
producing iteration-dependent thicknesses as given by equa- 
tion (11). From these constraint thicknesses, a set of active 
constraints is identified and a reduced set of “active” con- 
straints is determined. The reduced set of active constraints is 
obtained by considering the user-specified (1) minimum and 
maximum percentage of active constraints identified for each 
constraint type, (2) minimum and maximum number of 
active constraints to be used (by constraint type), and 
(3) overall minimums and maximums. Thus, fewer, or more, 
constraints than the number of active constraints identified 
from the constraint thickness criteria may be included in the 
reduced set for subsequent calculations. 
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For the reduced set of constraints, the user-specified pair 
of OC update methods (one for the Lagrange multipliers and 
another for the design variables) is invoked, using iteration- 
dependent acceleration parameters. Furthermore, if more than 
one pair of methods is chosen, a weighted average of the 
resulting updates is taken. This average is based on the 
weights also supplied by the user in the input data. The 
rescaling vector D is calculated after the updates to the 
Lagrange multipliers but prior to updates to the design 
variables. 

At the end of each iteration, the design is normally 
rescaled to ensure feasibility, unless otherwise specified by 
the user. This core iteration process continues either until 
convergence is achieved or the maximum number of itera- 
tions is reached. 



Figure 3. — Optimality criteria code levels and constraint inclusion 
strategies. (Strategy 0 allows the level to be skipped). 


The optimization procedure has three levels, and each 
level has several choices of strategies, as depicted in figure 3. 
The strategy at each level is specified by an input parameter. 

• Level I considers optimization for one constraint type, 
such as stress constraints (Strategy 1), displacement con- 
straints (Strategy II), or frequency constraints (Strategy III). 

• Level II allows simultaneous consideration of two types of 
constraints: displacement and frequency (Strategy I), stress 
and frequency (Strategy II), or stress and displacement 
(Strategy III). 

• Level III allows the simultaneous consideration of all 
three constraint types (i.e., stress, displacement, and 
frequency). 

Any levels may be skipped. For example, one could 
choose to solve a problem that contained all three types of 
constraints by skipping the first two levels and considering 
all constraints simultaneously at Level III. Alternatively, one 
could use Level I, Strategy I for stress constraints, then uti- 
lize these results as an initial design for Level II, Strategy III 
to include both stress and displacement constraints. Finally, 
these results could be used to initiate Level III, thus including 
the frequency constraints. This second alternative is generally 
recommended. 

Besides these three sets of iteration loops, the initial and 
final design are rescaled to ensure design feasibility. The final 
design may also be rescaled to ensure that at least one con- 
straint becomes active. In addition, options are available to 
specify convergence criteria, the maximum iteration limits, 
and the initial values of the Lagrange multipliers. Further- 
more, users can remove restrictions normally placed on the 
Lagrange multipliers and the design variables during each 
iteration by (1) allowing infeasibility of the design at most 
intermediate stages, (2) allowing the Lagrange multipliers to 
become negative, and (3) allowing the rescaling vector D to 
include negative values. 

C. Analysis Modules 

The behavior constraints and objective function are gener- 
ated in the analysis module. The objective function is given 
simply by the weight of the structure. The structure is ana- 
lyzed to calculate the behavior constraints following a finite 
element technique. A choice of three analysis methods are 
currently available: (1) the displacement (or stiffness) 
method, (2) the integrated force method, and (3) a simplified 
integrated force method. 

L Displacement method . — The stress, displacement, and 
frequency constraints for the structure under single or mul- 
tiple load conditions are obtained with the displacement 
(stiffness) method as implemented in the ANALYZ/DANLYZ 
code developed at Wright-Patterson Air Force Base (ref. 17). It 
generates displacement constraints at specified nodes and 
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directions and/or specified frequency constraints, along with 
all stress constraints. The design sensitivity matrix for stress, 
displacement, and frequency constraints is obtained analyti- 
cally. However, finite difference sensitivity calculations are 
also available as an option. 

The information passed on to the optimization module 
from the analysis module consists of the value of the objec- 
tive function, the gradient of the objective function, the con- 
straint array, and the gradient (or sensitivity matrix) of the 
constraints. Analysis and sensitivity calculations via the 
displacement (or stiffness) method have become routine pro- 
cedures, so no further discussion of this method is presented 
in this report (see refs. 18 and 19). 

It should be noted here that certain apparent anomalies are 
observed with certain constraint gradient information pro- 
vided by the ANALYZ/DANLYZ implementation of the 
analysis method. This is not deemed sufficient, however, to 
preclude the use of this code for comparing optimizers. 

2. Integrated force method. — The force method analysis 
generates and passes along the same basic information (the 
objective function, its gradient, the constraints, and their sen- 
sitivities) to the optimization module. Because the integrated 
force method is relatively new, a review of the basic equa- 
tions for the analysis and sensitivity calculations is presented 
(see ref. 20). 

The integrated force method (IFM), considers all the inter- 
nal forces F to be the simultaneous unknowns. The force 
equilibrium equations [B]F =P and the strain compatibility 
conditions [C][G]F =SR are concatenated to obtain the gov- 
erning equations of the method: 


The basic equation of IFM for frequency analysis of an 
undamped structure is given by 


[[S]-a> 2 [A/*]]/ m = 0 

where 


\m*] = 

[M j\g]] 

l J 

[0] J 


(28) 


(29) 


here, [A f] is the mass matrix, co is the circular frequency, and 
f m is termed the “force mode shape” of the eigenvalue 
problem. 

Note that in equation (26a), matrices [£] and [C] are inde- 
pendent of the design variables and need not be regenerated 
after the first analysis of a particular problem. Only the block 
diagonal flexibility matrix [G], which contains the design 
variables, has to be regenerated for reanalysis. 

The equations used for the sensitivity calculations are 
given for trusses. Here, the areas X are taken as the design 
variables. 

The sensitivity matrix, of dimension n X n> for stress con- 
straints is 


[Va] = [D] 7 '[ij-[-^-] (30) 


W 1; 

jp_~ 

Lmr 

~ SR 


or 

[S]F = F* (26b) 

where [F] is the m X n equilibrium matrix, [C] is the r X n 
compatibility matrix, [G] is the n X n concatenated flexibility 
matrix that links deformations to forces F as Q3 =[G]F ), 
P is the m component load vector, SR is the r component 
effective initial deformation vector (SR = ~[C]/?o)> Po IS 
the n component initial deformation vector, [S] is the n X n 
governing matrix, and m + r = n. The matrices [B], [C], [G], 
and [S] are banded and have full row ranks of m, r, n, and n , 
respectively. 

The solution of equation (26a) or (26b) yields the n forces 
F . The m displacements X are obtained from the forces as 

* = [7]([C]F + A)) (27) 


where 


[V<T] = [Va,,V<T 2 Va„] 


Vaf 


do t d(Tj 
dA x * dA 2 ’ 


dOj y 

) 


T 


[D] = [5]- 1 



and [1/A], [F/A 2 ], [G], and [F] are diagonal matrices of 
dimension n X n, whose diagonal elements are 1/A,-, Fj/A?, 
£i/(a} Ej), and F f , respectively. Here, ft and are the lengths 
and moduli of elasticity of each element, and F; is the ith 
component of the force vector F . 

The gradient matrix, of dimension n X m, for displace- 
ment constraints is 


Here, [J] is them X n deformation coefficient matrix defined T 

as ([J] = m rows of [ [S]" 1 ] 7 ). [VX] = \[J][G][D] + [J)[s dg ]\ (31) 
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The elements of the diagonal matrix [Sjg] are given by 



(31a) 


the illustration in appendix C, replacement of the expression 
[7I[G][D] in equation (31) by a null matrix, simplifies the 
displacement sensitivity, which becomes 

[vxr«toj/] r (33) 


where ga is the value of the iih diagonal element of the flex- 
ibility matrix [G]. 

The gradient matrix for frequency constraints may be 
obtained by following the approach of Rudisill (ref. 19 
and 20). A typical element of this sensitivity matrix is given 
by 


do 

dAi 


1 

2 co 



V 


-CD 2 


\ } 
f 

J 

fol J 

J 

L 



/ 


(32) 


where the superscript a represents an approximation to the 
sensitivity matrix. 

The stress gradient given by equation (30) can be approxi- 
mated by dropping the first term, which represents the effect 
of compatibility in the sensitivity calculations. The approxi- 
mate expression for the stress sensitivity has the following 
form: 



(34) 


The effects of the simplified sensitivity formulas are dis- 
cussed in section V. 


where 


D. Linking of Design Variables 



Design variables in a given problem can be reduced by 
linking them into groups. Specifically, the design variables 
are divided into a small number of groups. Linking factors 
can be assigned to elements of the group. 

Consider, for example, the ring structure depicted in fig- 
ure 4, The 60 element areas are divided into 25 linked groups 
(as depicted in table 23); 

Group 1 is {x49’XsOT-’X6o} = {l-QX\>l-Oxt 

Group 2 is {X\,Xn) = 1-0*2 } 

Group 3 is foau} = (l •0*}', 1 .0*3" } 

Group 25 is {* 36 ,* 48 } = {l.O^.l.O*^} 

where xf represents the yth linked design variable. The link- 
ing factors in this example are 1,0; these factors, however, 
can be different. 

The efficiency of this technique is tested in section IV and 
discussed in section V. 

E. CometBoards 


C { is the /th column of [C], fi and f r are the left and right 
eigenvectors of the eigenvalue problem, respectively, and 7; 
may be obtained through its relationship with [S], 

3 . Simplified integrated force method . — The sensitivity 
analysis equations of the integrated force method are simpli- 
fied, and alternative versions that are numerically less expen- 
sive than the original formulas are obtained. As suggested by 


The OC code was incorporated into the Comparative 
Evaluation Test Bed of Optimization and Analysis Routines 
for the Design of Structures (CometBoards). This test bed is 
available on the VM/CMS computer system, as well as on the 
Posix-based Cray and Convex computers and Iris worksta- 
tions. The command-level user interface is similar on all four 
platforms. The command syntax for invocation of the OC 
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Figure 4. — Sixty-bar trussed ring. (Elements are circled, nodes are not.) 


code on VM/CMS follows. Bold type should be input 
exactly as shown. Italic type highlights parameters where the 
user has a choice of options. 

optimize oc analyze r>fileid 


where 



analyzer is 

disp 

for the displacement method 
analyzer 


ifm 

for the integrated force method 
analyzer 


sifmsd 

for the simplified integrated force 
method analyzer 

file id is 

the file identification for the output file on 


VM CMS (the default is for output to come 
to the screen) 


The user on VM/CMS is prompted for the fileids of three 
input data files. For example, to run the optimality criteria 
method optimizer with the displacement analyzer on a user- 
specified problem with stress, displacement, and frequency 


constraints, and to place the output into the following File — 
problsdf output a — the user would type 

optimize oc disp>problsdf output a 

Then, following each of three prompts, the user would enter 
each of three fileids containing input data, such as 

problsdf anldat a 
problsdf idsdat a 
problsdf ocdat a 

The anldat file would include the specification of the finite 
element model. The idsdat file would contain design input 
data needed for constraint generation. The ocdat file would 
include parameters for the optimizer, including which 
optimality criteria update formulas are to be used. 

IV. Numerical Examples 

Several numerical examples for minimum weight design 
under single and multiple load conditions with stress, dis- 
placement and/or frequency constraints are solved using 
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the OC code. The results obtained are compared with a math- 
ematical programming technique referred to as a Sequence of 
Unconstrained Minimizations Technique (SUMT), which is 
available in the CometBoards test bed. Solutions obtained by 
SUMT are qualified by other optimization methods. Results 
from these example problems are briefly presented in this 
section. 

Seventeen independent combinations of update formulas 
within the OC code were used to solve the example prob- 
lems. These 17 procedures include optimality criteria meth- 
ods, modified OC methods, a fully utilized design, and 
hybrid methods, which combine some of the OC methods 
with a fully utilized design. The salient features of the 
17 methods are described next. An OC method, method 6 for 
example, updates the Lagrangian multipliers by using the 
exponential form, whereas it updates the design by using the 
reciprocal form. 

The first six methods are standard optimality criteria 
methods given in references 3 and 4. Each method combines 
one update formula for the calculation of the Lagrange mul- 
tipliers and one formula for updates to the design variables. 
These combinations are given in the following chart. 


The next seven methods are modified OC methods. Their 
combinations of update formulas are given in the following 
chart. 


Method 

Lagrange multiplier 
update formula 

Design variable 
update formula 

7 

Diagonalized Lagrange inverse 

Exponential form 

8 

Diagonalized Lagrange inverse 

Linearized form 

9 

Diagonalized Lagrange inverse 

Reciprocal form 

10 

Unrestricted Lagrange inverse 

Exponential form 

11 

Unrestricted Lagrange inverse 

Linearized form 

12 

Unrestricted Lagrange inverse 

Reciprocal form 

13 

Lagrange inverse form 

Melange form 


The next three methods are hybrid methods, each using a 
hybrid design variable update formula as described in 
section II. The following chart depicts the update formula 
combinations for each of these methods. 


Method 

Lagrange multiplier 

Design variable 


update formula 

update formula 

14 

Lagrange inverse form 

Hybrid exponential form 

15 

Lagrange inverse form 

Hybrid linearized form 

16 

Lagrange inverse form 

Hybrid reciprocal form 


Method 17 is the fully utilized design for stress, displace- 
ment, and frequency constraints. 

Seven sets of examples consisting of a total of 31 prob- 
lems are given. The problems range from a modest truss 
example with 5 design variables under displacement con- 
straints to a difficult intermediate complexity wing problem 
with different element types, 320 constraints, and 57 linked 
design variables. The load conditions, mass distributions, 
stress, displacement, and frequency constraints are chosen to 
ensure that, at optimum, several behavior constraints are 
active. For example, one three-bar truss problem results in a 
total of nine active constraints at optimum, consisting of one 
frequency, two displacement, and six stress constraints. In all 
cases, the objective is to minimize the weight of the structure. 

A. Five-Bar Truss 

The five-bar truss reported in reference 4 consists of a set 
of four problems. The truss is made of aluminum, with a 
Young’s modulus Eof 10 000 ksi, a Poisson’s ratio vof 0.3, 
and a weight density p of 0.1 lb/in. 3 Two related problems 
are also included in this first example set. Displacement limi- 
tations are the only behavior constraints. The structures are 
depicted in figures 5 and 6, whereas the loads and constraints 
are given in tables 1 and 2, respectively. Note that the 



Figure 5. — Five-bar truss. (Elements are circled, nodes are not.) 


Method 

Lagrange multiplier 
update formula 

Design variable 
update formula 

1 

Linear form 

Exponential form 

2 

Linear form 

Linearized form 

3 

Linear form 

Reciprocal form 

4 

Exponential form 

Exponential form 

5 

Exponential form 

Linearized form 

6 

Exponential form 

Reciprocal form 
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Figure 6.— Modified five-bar truss. (Elements are circled, nodes 
are not.) 


boundary conditions differ between the original four prob- 
lems and the last two. The examples are designated as 
BarSL.l, Bar5L.2, Bar5L.3, Bar5L.4, Bar5L.2a, and 
Bar5L.2b. The optimal solutions given by SUMT show that 
all (either one or two) displacement constraints become 
active for all the problems except problem Bar5L.2, where 
only one of the two becomes active. 

Solutions for the six problems are obtained using all 
17 methods (OC-type, hybrid, and FTJD), as well as SUMT, 
and are presented in tables 3 to 8. The results for five of the 
OC methods (4, 5, 12, 13, and 16) agree well with the SUMT 


TABLE I.— FIVE-BAR TRUSS 
EXAMPLES: LOAD SPECIFICATIONS 


Problem 

Number of 
boundary 
conditions 

Load components, 
kips 

Node 

¥ 

X 

¥ 

y 

Bar5L. I 


X 

4 

0 

0 

100.0 

Bar5L.2 



4 



100.0 

Bar5L.3 



2 



-100.0 




4 



1 

30.0 

Bar5L.4 



2 







¥ 

4 


? : 

1 

f 

Bar5L.2a 


3 

2 

0.0 

100.0 

Bar5L.2b 


3 

2 

0.0 

100.0 


TABLE 2. - FIVE-BAR TRUSS EXAMPLES: 
CONSTRAINT SPECIFICATIONS 


Problem 

Constrain! 

type 

Description of 
constraints along 
y direction 

Node 

Magnitude, 

in. 

Bar5L. 1 


4 

2.0 

Bar5L.2 


2 

6.0 



4 

6.0 

Bar5L.3 

Displacement 

2 

-.1 



4 

.1 

Bar5L.4 


2 

2.0 



4 

2.0 

Bar5L.2a 


4 

2.0 

Bar5L.2b 

Displacement 

2 

6.0 



4 

2.0 

Bar5L. 1 




Bar5L.2 




Bar5L.3 

Minimum 

Xi = 

: 0.001 in. 2 

Bar5L.4 

area 

0 = 

: 1,2,... ,5) 

Bar5L.2a 




Bar5L.2b 





TABLE 3.- FIVE-BAR TRUSS EXAMPLES: 
PROBLEM BAR5L.1 
[Displacement constraints only.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 2 

X 

X 

Xi 

Xi 

X3 

X4 

Xs 

1 

1 

1 

45.016 

0.001 

1.499 

0.001 

2.120 

0.001 

2 

1 

2 

45.016 

.001 

1.499 

.001 

2.120 

.001 

3 

1 

3 

79.970 

1.998 

.001 

2.826 

.003 

1.998 

4 

2 

1 

45.016 

.001 

1.499 

.001 

2.120 

.001 

5 

2 

2 

45.016 

.001 

1.499 

.001 

2.120 

.001 

6 

2 

3 

79.970 

1.998 

.001 

2.826 

.003 

1.998 

7 

3 

1 

45.034 

.001 

1.500 

.001 

2.121 

.001 

8 

3 

2 

45.034 

.001 

1.500 

.001 

2.121 

.001 

9 

3 

3 

45.034 

.001 

1.500 

.001 

2.121 

.001 

10 

4 

1 

45.040 

.001 

1.524 

.001 

2.105 

XX>1 

11 

4 

2 

45.171 

.001 

1.624 

.001 

2.043 

.001 

12 

4 

3 

45.171 

.001 

1.624 

.001 

2.043 

.001 

13 

5 

4 

45.268 

.015 

1.488 

.022 

2.105 

.015 

14 

5 

5 

45.040 

JCX)1 

1.524 

.001 

2.105 

.001 

15 

5 

6 

45.171 

.001 

1.624 

.001 

2.043 

.001 

16 

5 

7 

45.171 

.001 

1.624 

.001 

2.043 

.001 

17 

0 

8 

62.228 

1.068 

1.068 

1.068 

1.060 

1.068 

SUMT" 

45.029 

.001 

1.501 

.001 

2.119 

.001 


"Sequence of Unconstrained Minimizations Technique. 
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TABLE 4. - FIVE-BAR TRUSS EXAMPLES: 
PROBLEM BAR5L2 
[Displacement constraints only.] 


Methods 

Optimality 

criteria 

(OQ 

formulas 

Optimum 

weight, 

H> 

Optimum design variables, 
in. 

X 

X 

Xi 

X2 

X3 

X4 

X5 

1 

1 

1 

20.743 

0356 

0356 

0356 

0356 

0356 

2 

1 

2 

20.743 

.356 

356 

356 

356 

356 

3 

1 

3 

20.743 

356 

356 

356 

356 

.356 

4 

2 

1 

15.016 

.001 

.499 

.001 

.706 

.001 

5 

2 

2 

15.016 

.001 

.499 

.001 

.706 

.001 

6 

2 

3 

15.016 

.001 

.499 

.001 

.706 

.001 

7 

3 

1 

15.034 

.001 

300 

.001 

.707 

.001 

8 

3 

2 

15.034 

.001 

300 

.001 

.707 

.001 

9 

3 

3 

15.034 

.001 

300 

.001 

.707 

.001 

10 

4 

1 

15.036 

.001 

.508 

.001 

.702 

.001 

11 

4 

2 

15.080 

.001 

341 

.001 

.681 

.001 

12 

4 

3 

15.080 

.001 

.541 

.001 ! 

.681 

.001 

13 

5 

4 

15.084 

.005 

.496 

.007 

.702 

.005 

14 

5 

5 

15.036 

.001 

308 

.001 

.702 

.001 

15 

5 

6 

15.080 

.001 

341 

.001 

.681 

.001 

16 

5 

7 

15.080 

.001 

341 

.001 

.681 

.001 

17 

0 

8 

20.743 

356 

356 

356 

356 

356 

SUMT* 

15.026 

.001 

302 

.002 

.704 

.001 


‘Sequence of Unconstrained Minimizations Technique. 


TABLE 6. -FIVE-BAR TRUSS EXAMPLES: 
PROBLEM BAR5L.4 
[Displacement constraints only.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 

X 

X 


Xi 

12 

X.i 

X4 

Xs 

1 

1 

, 

90.010 

1.500 

1500 

2.121 

2.121 

0.001 

2 

1 

2 

90.010 





2.121 

2.121 

.001 

3 

1 

3 

90.087 





2.122 

2.122 

.009 

4 

2 

1 

90.010 





2.121 

2.121 

.001 

5 

2 

2 

90.010 





2.121 

2.121 

.001 

6 

2 

3 

90.023 



1 


2.121 

2.121 

.002 

•7 

3 

1 



_ 



- 

— 

_ 

-- 

- 

-- 

- 

— 

8 

3 

2 

97.982 

1.485 

1.769 

2.100 

2.502 

.034 

9 

3 

3 

96.051 

1.702 

1.500 

2.407 

2. 

21 

. 

301 

10 

4 

1 

90.010 

1500 



2.121 





11 

4 

2 

90.010 











12 

4 

3 

90.010 











13 

5 

4 

90.060 









.006 

14 

5 

5 

90.010 









.001 

15 

5 

6 

90.010 









B01 

16 

5 

7 

90.010 









.001 

17 

0 

8 

111368 

1.914 

1.914 

1.914 

1.914 

1.914 

SUMT* 

90.063 

1.498 

1.498 

2.124 

2.124 

.002 


‘Method failed for this problem. 

^Sequence of Unconstrained Minimizations Technique. 


TABLE 5. - FIVE-BAR TRUSS EXAMPLES: 
PROBLEM BAR5L3 
[Displacement constraints only.] 


Methods 

Optimality 

criteria 

(OQ 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in* 

X 

X 

Xi 

X2 

X3 

X4 

Xs 

1 

l 

1 

50.043 

0.001 

0.001 

0.001 

0.001 

5.000 

2 

1 

2 

50.043 

.001 

.001 

.001 

.001 

5.000 

*3 

1 

3 

— 

— 

— 

— 

— 

- 

-- 

4 

2 

1 

50.047 

.001 

.001 

.001 

.001 

5.000 

5 

2 

2 

50.047 

.001 

.001 

.001 

.001 

5.000 

*6 

2 

3 

— 

— 

— 

— 

— 

- 

-- 

7 

3 

1 

53.065 

.065 

.065 

.065 

.065 

4.991 

8 

3 

2 

55.971 

.122 

.122 

.131 

.131 

4.983 

9 

3 

3 

60.421 

.216 

216 

.226 

.226 

4.971 

10 

4 

1 

50.000 

BOO 

.000 

.000 

.000 

5.000 

11 

4 

2 

50.000 

.000 

BOO 

.000 

.000 



12 

4 

3 

50.000 

.000 

.000 

.000 

.000 

1 



13 

5 

4 

50.073 

.001 

.001 

.002 

.002 



14 

5 

5 

50.000 

.000 

.000 

.000 

.000 



15 

5 

6 

50.000 

BOO 

.000 

.000 

.000 



16 

5 

7 

50.000 

.000 

.000 

.000 

.000 



17 

0 

8 

257.758 

4.422 

4.422 

4.422 

4.422 

4.422 

SUMT b 

50560 

.008 

.008 

.013 

.013 

5.002 


‘Method failed for this problem, 

•^Sequence of Unconstrained Minimizations Technique. 


TABLE 7. -FIVE-BAR TRUSS EXAMPLES: 
PROBLEM BAR5L.2a 
[Displacement constraints only.] 


Methods 

Optimality 

criteria 

(OQ 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 

X 

X 


Xi 

12 

X3 

X4 

Xj 

1 

1 

1 

45.034 

0.001 

1.500 

0.001 

2.121 

0.001 

2 

1 

2 

45.034 

.001 

1.500 

.001 



.001 

3 

1 

3 

45.106 

.003 

1.499 

.003 



.003 

4 

2 

1 

45.034 

.001 

1.500 

.001 



.001 

5 

2 

2 

45.034 











6 

2 

3 

45B42 











7 

3 

1 

45.036 











8 

3 

2 

45.035 











9 

3 

3 

45.035 











10 

4 

1 

45.036 











11 

4 

2 

45B35 











12 

4 

3 

45.035 






f 





13 

5 

4 

45.078 

.002 



.002 



J 

002 

14 

5 

5 

45.036 

.001 



.001 



.001 

15 

5 

6 

45B35 

.001 



.001 



.001 

16 

5 

7 

45.035 

.001 



.001 



.001 

17 

0 

8 

111368 

1.914 

1.914 

1.914 

1.914 

1.914 

SUMT 

45B51 

.001 

1.497 

0.001 

Z124 

.001 


‘Sequence of Unconstrained Minimizations Technique. 
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TABLE 8. - FIVE-BAR TRUSS EXAMPLES: 
PROBLEM BAR5L.2b 
[Displacement constraints only.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 2 



\ 

X 

Xi 

X2 

X3 

X4 

Xs 

1 

1 

1 

111.568 

1.914 

1.914 

1.914 

1.914 

1.914 

2 

1 

2 

103.459 

.197 

3.223 

.197 

4.600 

.197 

3 

1 

3 

67.518 

.225 

1532 

.225 

3.148 

.225 

4 

2 

1 

47540 

.001 

1.500 

.001 

2.121 

.252 

5 

2 

2 

47.539 

.001 

1500 

.001 

2.121 

.251 

"6 

2 

3 

— 

— 

— 

— 

— 

— 

*7 

3 

1 

— 





— 





8 

3 

2 

53.208 

.001 

1.694 

.001 

2.395 

.236 

9 

3 

3 

48.446 

.001 

1.500 

.001 

2.121 

.341 

10 

4 

1 

60.217 

.009 

1500 

.009 

2.121 

1.500 

11 

4 

2 

76.593 

.002 

2.482 

.002 

3510 

.209 

12 

4 

3 

47.551 

.001 

1.500 

.001 

2.121 

.253 

13 

5 

4 

48.368 

.001 

1.528 

.001 

2.162 

.248 

14 

5 

5 

60.217 

.009 

1.500 

.009 

2.121 

1500 

15 

5 

6 

76593 

.002 

2.482 

.002 

3.510 

.209 

16 

5 

7 

47551 

.001 

1500 

.001 

2.121 

.253 

17 

0 

8 

111568 

1.914 

1.914 

1.914 

1.914 

1.914 

SUMT* 

47540 

.001 

1.496 

.001 

2.125 

.250 


"Method failed for this problem. 

b Sequence of Unconstrained Minimizations Technique. 


answers for all six problems (although method 13 shows a 
1.7 percent difference for the last problem). Four additional 
methods (10, 11, 14, and 15) show results that agree with 
SUMT for all but the last problem (Bar5L.2b). Among most 
of the other OC and hybrid methods, agreement is seen for at 
least half the problems. The fully utilized design approach is 
not applicable for this problem because no stress constraints 
are specified, but the results are included for completeness. 

B. Tapered Five-Bar Truss Under Stress, Displacement, 
and Frequency Constraints 

A tapered five-bar truss, shown in figure 7, provides three 
more example problems. It is made of steel, with a Young’s 
modulus E of 30 000 ksi, a Poisson’s ratio v of 0.3, and a 
weight density p of 0.284 lb/in. 3 The truss is subjected to two 
load conditions and has five stress and two displacement con- 
straints per load condition, along with a frequency constraint, 
as shown in tables 9 and 10. These minimum weight design 
problems are designated by (1) Bar5.a, where only stress con- 
straints are specified (of which eight are active with SUMT’s 
optimal solution) and the results are depicted in table 11; 
(2) BarS.b, where stress and displacement constraints are 
specified (of which three displacement and six stress con- 
straints are active at optimum) and the results are shown in 
table 12; and (3) Bar5.c, where stress, displacement, and fre- 
quency constraints are all specified (of which one frequency, 
two displacement, and two stress constraints are active at the 
optimal solution) and the results are given in table 13. 


fy 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

100 |- 



Figure 7.— Tapered five-bar truss. (Elements are circled, nodes are not.) 
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TABLE 9. -TAPERED FIVE-BAR TRUSS 
PROBLEM: LOAD SPECIFICATIONS 


Problem 

Load 

conditions 

Load components, 
kips 



Node 

p* 

p y 


1 

2 

120.0 

120.0 

Bar5.a 

4 

35.0 

25.0 

BarS.b 

BarS.c 

11 

2 

0.0 

-65.0 


4 

0.0 

-75.0 

Bar5.c 

Lumped 

masses 

2 

4 

m x = m y = 114 lb 

m x - m y ~ 165 lb 


TABLE 10. -TAPERED FIVE- BAR TRUSS: 
CONSTRAINT SPECIFICATIONS 


Problem 

Constraint 

type 

Constraint description 

Bar5.a 

Stress 

Minimum area 

a,£a 0 (* = 1,2,... ,5) 

o 0 = 20 ksi 

x, ^ 0.25 in. 2 (i= 1,2,...,5) 

Bar5.b 

Stress 

Displacement 
Minimum area 

Same as Bar5.a 
Constraints along y direction 
1.9 in. at node 2 
1.9 in. at node 4 
Same as BarS.a 

BarS.IIc 

Stress 

Displacement 

Frequency 

Same as Bar5.a 
Same as Bar 5.b 
fzfolfi b=15Hz 


Minimum area 

Xi 2 4.0 in. 2 (r = 1,2,.. .4) 


TABLE 11. -TAPERED FIVE-BAR TRUSS: 
PROBLEM BAR5.a 
[Stress constraints only.] 



t* 
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Figure 8.— Ten-bar truss. (Elements are circled, nodes are not.) 
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From tables 11 to 13, it is observed that method 17, the 
fully utilized design, yields solutions close to the SUMT 
results for all the three problems. Methods 5, 6, and 13 show 
satisfactory performance with less than a 3 percent difference 
with SUMT (except for method 6, which differs by 5.6 per- 
cent in the first problem, Bar5.a). Performance of the three 
hybrid methods differs depending on the type of constraints 
considered. For stress constraints only (Bar5.2a), results 
obtained by all three hybrid methods are good, differences in 
the optimal weight being less than 1 percent; but results are 
poor when displacement and/or frequency constraints are 
included. The performance of the other 10 OC methods for 
this set of problems is generally poor. 

C. 10-Bar Truss Design 

The 10-bar aluminum truss, shown in figure 8, is taken as 
the third example set. In the first two cases, the truss is sub- 
jected to two load conditions, with 10 stress and 2 displace- 
ment constraints per load condition, along with a frequency 
constraint, (for a total of 25 constraints), as given in tables 14 
and 15. Two design situations are considered. Problem 
BarlO.a has all 10 bar areas considered as independent design 
variables. The second case, BarlO.b, is obtained by linking 
the variables into a set of five independent linked design 


TABLE 12. -TAPERED FIVE-BAR TRUSS: 
PROBLEM BAR5.b 
[Stress and displacement constraints.! 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 

X 

X 

Xi 

X2 

X3 

X4 

Xi 

1 

1 

1 

7 743.920 

20.001 

20.001 

20.001 

20.001 

20.001 

2 

1 

2 

6 997.886 

20365 

17.657 

17.073 

18.400 

10.433 

3 

1 

3 

6 985.157 

20.420 

17580 

16.817 

17.999 

14.018 

4 

2 

1 

9 266.655 

43.652 

7.979 

41588 

5.200 

5.200 

5 

2 

2 

6 632.620 

23391 

13.068 

18.703 

15.133 

5557 

6 

2 

3 

6 763.739 

21.412 

16332 

15.782 

18.232 

5.489 

7 

3 

1 

12 870.526 

5.999 

63.995 

.080 

66.317 

13510 

8 

3 

2 

8 561.541 

22.726 

28.900 

9588 

30.009 

5.487 

•9 

3 

3 


— 

— 

— 

— 

— 

10 

4 

1 

7 456582 

31.803 

8.796 

29.071 

9,894 

2.053 

11 

4 

2 

7 145541 

29.424 

9507 

26.777 

10568 

1.852 

M2 

4 

3 


— 

— 

i 

— 

— 

13 

5 

4 

6 549.017 

26.469 

10330 

21329 

11.973 

1.456 

14 

5 

5 

14 873.955 

67.303 

13.761 

61.949 

15.913 

2377 

15 

5 

6 

17 976.304 

82.470 

15521 

76.026 

18.108 

2.453 

M6 

5 

7 


— 

— 

— 


— 

17 

0 

8 

6 549.674 

27376 

9.421 

21.909 

11.266 

1.633 

SUMT* 

6 541559 

26.208 

10397 

21327 

11.931 

1.809 


■Method failed for this problem. 

"Sequence of Unconstrained Minimizations Technique. 


TABLE 13. -TAPERED FIVE-BAR TRUSS: 
PROBLEM BAR5.C 

[Stress, displacement, and frequency constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 

in* 

X 

X 

Xl 

X2 

X3 

Xa 

X5 

1 

1 

1 

8083.500 

20.878 

20.878 

20.878 

20.878 

20.878 

2 

1 

2 

7 391.503 

20572 

18.983 

18.753 

19.419 

10351 

3 

1 

3 

7 377393 

21.088 

18.678 

18555 

19.044 

14341 

4 

2 

1 

8083500 

20.878 

20.878 

20378 

20378 

20378 

5 

2 

2 

6890.76! 

23.498 

14.643 

17.974 

16578 

8.143 

6 

2 

3 

6892336 

21.708 

16.493 

16.206 

18311 

6.111 

*7 

3 

1 




— 

— 


— 

■8 

3 

2 


— 

— 

— 

— 

— 

■9 

3 

3 


— 

— 

— 



— 

MO 

4 

1 




— 

— 

— 

— 

Ml 

4 

2 


— 

— 

— 

— 

— 

12 

4 

3 

8 521.432 

14.827 

32.046 

9.998 

33.685 

6.590 

13 

5 

4 

6798.956 

28586 

8.747 

24.184 

10.426 

4.126 

14 

5 

5 

13 512.923 

57520 

16.029 

52576 

18.021 

4.000 

15 

5 

6 

18037.263 

81.186 

17.023 

74.704 

19.638 

4.000 

M6 

5 

7 


— 

— 

— 


— 

17 

0 

8 

6 742305 

26.898 

10.664 

21511 

12591 

4.141 

SUMT" 

6 958.113 

22.492 

15.848 

19.740 

15.937 

4.000 


■Method failed for this problem. 

"Sequence of Unconstrained Minimizations Technique. 


TABLE 14. -BAR TRUSS: LOAD 
SPECIFICATIONS 


Problem 

Load 

conditions 

Load components, 
kips 


D 

wm 





120.0 





60.0 


I 


175 

12.5 



5 

175 

25.0 




17.5 

25.0 

BarlO.a 





BarlO.b 




-50.0 





-25.0 


II 



-375 



5 


-75.0 




0.6 

6.0 

BarlO.c 



6.0 

6.0 

BarlO.d 



1.75 

1.25 



5 

.17 

2.5 

BarlO.a 


B 


K * 751b 

BarlO.b 

Lumped 

n 

BBS 

v m 135 lb 

BarlO.c 

masses 

mm 

BBS 

75 lb 

BarlO.d 


B 


y = 135 lb 
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TABLE 15.- 10-BAR TRUSS: 
CONSTRAINT SPECIFICATIONS 


TABLE 16. -10-BAR TRUSS: DESIGN 
VARIABLE LINKAGE 


Problem 

Serial 

Design 

Members 


number 

variable 

linked 


n 


1,5 



1 

2,4 

BarlO.b 

H 

jrcH 

3,6 




7,8 



B 

9,10 


1 

1 

1 



■ I 

2,4 

BarlO.d 


1 

5 




3,9,10 


5 

m 

6,7,8 


Problem 

Constraint 

type 

Constraint description 


Stress 

C^Oq (i= 1,2,... ,10) 
a 0 - 10 ksi 

BarlO.a 
Bar lO.b 

Displacement 

Constraints along y direction of magnitude: 
2.2 in. at node 3 
2.2 in. at node 4 


Frequency 

/;>/□;/□ = 26 Hz 


Minimum area 

Xi = 4.0 in. 2 (t = 1,2,..., 10) for BarlO.a 

Xi = 4.0 in. 2 (i = 1,2,. ..,5) for BarlO.b 


Displacement 

Constraints along y direction of magnitude: 
1.5 in. at node 3 
1.5 in. at node 4 

BarlO.c 

BarlO.d 

Frequency 

/2>/ 0 ;/ 0 =15Hz 


Minimum area 

X,- = 0. 1 0 in. 2 (i'=l,2 10) for BarlO.c 

Xi = 0.10in. 2 (i = 1,2,.. .,5) for BarlO.d 


TABLE 17. — OPTIMUM DESIGN OF A 10-BAR TRUSS: PROBLEM BARlO.a 
[Stress, displacement, and frequency constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Normalized 

weight 

Normalized design variables 

X 

X 

Xi 

X2 

X3 

X4 

X5 

X6 

X7 

Xh 

X9 

Xio 

1 

1 

1 

1.78 

1.27 

1.69 

7.50 

4.43 

1.69 

7.50 

1.21 

1.47 

2.32 

3.82 

2 

1 

2 

1.94 

1.96 

.79 

3.35 

2.03 

2.43 

3.37 

1.97 

2.12 

1.06 

1.77 

3 

1 

3 

1.25 

1.08 

.89 

3.90 

2.32 

1.27 

3.93 

1.12 

1.15 

1.22 

2.01 

4 

2 

1 

1.23 

.89 

.79 

3.51 

2.07 

1.64 

3.51 

.57 

1.71 

1.09 

1.79 

5 

2 

2 

1.11 

1,04 

.69 

1.96 

1.96 

1.23 

2.32 

1.04 

1.14 

.61 

1.91 

6 

2 

3 

1.07 

1.05 

,81 

1.49 

1.55 

1.11 

1.49 

1.09 

1.03 

.77 

1.61 

7 

3 

1 

3.66 

3.38 

4.81 

1.00 

3.50 

3.48 

LOO 

3.63 

2.89 

5.89 

3.80 

8 

3 

2 

3.69 

3.41 

4.83 

1.00 

3.53 

3.48 

LOO 

3.77 

2.87 

5.92 

3.83 

•9 

3 

3 

Failed 

(a) 

(a) 

(a) 

W 

(a) 

(a) 

(a) 

— 

— 

— 

10 

4 

1 

14.06 

13.51 

19.32 

4.87 

13.11 

13.60 

5.80 

11.87 

10.73 

23.65 

14.20 

11 

4 

2 

2.28 

2.09 

3.19 

1.00 

1.84 

2.26 

LOO 

1.92 

1.91 

3.92 

2.04 

M2 

4 

3 

Failed 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

13 

5 

4 

1.18 

.89 

1.65 

1.08 

.64 

1.35 

1.82 

.73 

1.51 

2.06 

.55 

14 

5 

5 

1.70 

1.56 

2.29 

1.00 

1.19 

1.85 

1.15 

1.14 

1.62 

2.86 

1.32 

15 

5 

6 

1.69 

1.55 

2.29 

1.00 

1.18 

1.84 

1.15 

1.14 

1.62 

2.85 

1.32 

M6 

5 

7 

Failed 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

17 

0 

8 

1.05 

1.15 

1.18 

1.11 

.79 

.80 

1.21 

1.35 

.87 

1.18 

.95 

SUMT* 1 


1.00 

1.00 

1.00 

LOO 

1.00 

LOO 

LOO 

LOO 

1.00 

LOO 

1.00 

Weight, lb 


3 326.93 

56.53 

17.70 

4.00 

6.77 

46.10 

4.00 

24.72 

23.76 

12.91 

| 7.85 


a Metbod failed for this problem, 

b Sequence of Unconstrained Minimizations Technique. 
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groups as shown in table 16. At optimum, for both these 
problems, the SUMT results show that one frequency, two 
displacement, and seven stress (six, in the BarlO.b case) con- 
straints are active. 

Both. cases are solved by all 17 methods and SUMT. The 
results obtained are normalized with respect to the SUMT 
answers, and are given in tables 17 and 18. Overall, all the 
methods perform at about the same level for both problems 
(BarlO.a and BarlO.b). The performance of these methods for 
the 10-bar truss is similar to that of the last two 5-bar truss 
problems (Bar5.b and Bar5.c). Methods 5 and 17 yield results 
that differ by around 10 percent or less, compared with 
SUMT. Method 6 shows a 7 percent difference with problem 
BarlO.a, but 17 percent with BarlO.b, whereas method 13 
shows an 18 percent difference with BarlO.a, but less than 
1 percent difference with BarlO.b. 

One effect of linking of design variables (see 
section ELD) is illustrated by also subjecting this 10-bar truss 


TABLE 18.— 10-BAR TRUSS: PROBLEM BARI O.b LINKED 
DESIGN CONFIGURATION 
[Stress, displacement, and frequency constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Normalized 

weight 

Normalized design variables 

X 

X 

Xi 

X2 

X3 

X4 

X5 

1 

1 

1 

1.69 

1.42 

2.12 

5.34 

1.29 

2.69 

2 

1 

2 

1.31 

1.21 

1.39 

3.29 

1.14 

1.75 

3 

1 

3 

1.60 

1.29 

2.02 

5.02 

1.31 

2.55 

4 

2 

1 

1.83 

1.32 

2.56 

6.44 

1.47 

3.25 

5 

2 

2 

1.08 

1.09 

.99 

1.40 

1.08 

1.09 

6 

2 

3 

1.17 

1.16 

.97 

1.58 

1.21 

1.31 

a 7 

3 

1 










a 8 

3 

2 

— 

— 



— 


39 

3 

3 

— 

— 



— 


a 10 

4 

1 









fia 

a ll 

4 

2 

— 

— 

HB 


— 

BB 

a 12 

4 

3 

— 

— 

— 

— 

— 

— 

13 

5 

4 

1.00 

.98 

1.01 

1.02 

1.02 

1.05 

14 

5 

5 

1.55 

1.54 

1.80 

.93 

1.24 

2.17 

15 

5 

6 

2.18 

2.16 

2.52 

1.20 

1.74 

3.04 

a 16 

5 

7 

— 

— 

— 

— 

— 

— 

17 

0 

8 

1.00 

.98 

1.01 

1.02 

1.02 

1.01 

SUMT* 


1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

Weight, lb 


3 572.53 

53.92 

14.37 

5.72 

25.37 

11.33 


a Method failed for this problem. 

Sequence of Unconstrained Minimizations Technique. 


to a single load condition (see table 14), with displacement 
and frequency constraints (see table 15). As before, two 
design situations are considered. Problem BarlO.c has all 
10 bar areas considered as independent design variables. The 
last case, BarlO.d, is obtained by linking the variables into a 
set of five independent linked design groups as shown in 
table 16. At optimum, for these last two problems, the SUMT 
results show that one frequency and two displacement con- 
straints are active. 

Once again, in addition to SUMT, all 17 methods are used 
to solve these problems. The results are given in tables 19 
and 20. Methods 1 to 6 perform well, showing less than a 
1 percent difference with SUMT for both problems. Method 
13 shows less than a 2 percent difference. Methods 7, 10, 11, 
14, and 15 also perform well for Problem BarlO.c but do 
poorly for the linked case (BarlO.d). Additional discussion 
regarding the effect of linking of design variables is provided 
in section V. 

D. 60-Bar Trussed Ring 

A 60-bar trussed ring made of aluminum (fig. 4) is sub- 
jected to three load conditions (table 21). The constraints 
considered are specified in table 22. The 60 element areas of 
the structure are linked into 25 groups each considered as a 
design variable (table 23). The problem is solved for the fol- 
lowing six situations: (a) Bar60.a, stress constraints only; 
(b) Bar60.b, displacement constraints only; (c) Bar60.c, fre- 
quency constraints only; (d) Bar60.d, both stress and dis- 
placement constraints; (e) both displacement and frequency 
constraints; and (0 all three types of constraints. The chart 
below shows the number of active constraints of each type in 
SUMT’s solutions for all six cases. 


Number of active constraints by type 

Problem 

Stress 

Displacement 

Frequency 


30 

__ 

.. 


... 

1 

- 


... 

- 

1 


24 

1 

- 


— 

1 

1 

f 

17 

1 

1 


Results obtained for the six problems by all 17 methods, 
normalized by the SUMT results are given in table 24. For 
problem Bar60.b (displacement constraints only), almost all 
of the OC-type methods (1 to 16) perform well, with differ- 
ences from SUMT not exceeding 5 percent. For problem 
Bar60.c (frequency constraints only), 12 of the 16 OC-type 
methods perform well with differences of about 5 percent. 
Methods 4, 5, 6, and 13 exhibit larger differences. For prob- 
lem Bar60.e (displacement and frequency constraints), 
11 methods show less than 10 percent differences with 
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TABLE 19. — OPTIMUM DESIGN OF A 10-BAR TRUSS: PROBLEM BARIO.c 
[Displacement and frequency constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 2 

X 

X 

Xi 

X2 

X3 

X4 

Xs 

Xfi 

X? 

Xs 

X9 

Xto 

1 

1 

1 

414.064 

5.609 

3.265 

0.144 

0.514 

7.114 

0.144 

1.227 

3.282 

3.127 

0.611 

2 

1 

2 

413.962 

5.616 

3.260 

.137 

.517 

7.109 

.137 

1.232 

3278 

3.127 

.607 

3 

1 

3 

413.702 

5.630 

3.210 

.131 

.542 

7.116 

.131 

1.249 

3272 

3.086 

.629 

4 

2 

1 

412.263 

5.943 

2.388 

.103 

1.014 

7.178 

.103 

1.620 

3.097 

2.382 

1.048 

5 

2 

2 

412.242 

5.884 

2.498 

.106 

.943 

7.188 

.106 

1.553 

3.140 

2.480 

.981 

6 

2 

3 

412.186 

5.878 

2.500 

.106 

.939 

7.193 

.106 

1.548 

3.145 

2.483 

.978 

7 

3 

1 

415.584 

5.511 

3.018 

.353 

.285 

7.817 

.134 

1.020 

3.832 

2.913 

.390 

8 

3 

2 

506.694 

7.210 

2.952 

.100 

.469 

9254 

.100 

1.747 

4.300 

3.458 

.920 

9 

3 

3 

586.999 

8.933 

3.163 

.100 

1.061 

10.668 

.100 

2.579 

4.805 

2.875 

1.076 

10 

4 

1 

413.242 

5.948 

2.454 

.101 

.983 

7.172 

.101 

1.610 

3.105 

2.436 

1.022 

11 

4 

2 

413.238 

5.950 

2.426 

.101 

.984 

7.170 

.101 

1.612 

3.103 

2.434 

1.024 

M2 

4 

3 

Failed 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

13 

5 

4 

412.942 

6.023 

2.388 

.101 

1.045 

7.088 

.101 

1.696 

3.019 

2.366 

1.086 

14 

5 

5 

413.242 : 

5.948 

2.454 

.101 

.983 

7.172 

.101 

1.610 

3.105 

2.436 

1.022 

15 

5 

6 

413.238 

5.950 

2.453 

.101 

.984 

7.170 

.101 

1.612 

3.103 

2.434 

1.024 

M6 

5 

7 

Failed 

— 

— 

— 

— 

— 

— 

— 

— 

— 

— 

17 

0 

8 

771.019 

5.279 

5.279 

5.279 

5.279 

5279 

5.279 

5.279 

5.279 

5279 

5.279 

SUMT* 

411.796 

6.228 

2.009 

.103 

1.287 

7.005 

.104 

1.910 

2.872 

2.035 

1295 


•Method failed for this problem. 

b Sequence of Unconstrained Minimizations Technique. 
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TABLE 20. -OPTIMUM DESIGN OF A 10-BAR TRUSS: PROBLEM BARlO.d 
LINKED DESIGN CONFIGURATION 
[Displacement and frequency constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design variables, 
in. 2 

X 

X 

Xi 

X2 

X3 

X4 

X5 

Xft 

X? 

X8 

X9 

Xio 

1 

1 

1 

444.351 

6.851 

1.945 

1.770 

1.945 

6.595 

2.224 

2.224 

2.224 

1.770 

1.770 

2 

1 

2 

444.260 

6.853 

1.939 

1.766 

1.939 

6.603 

2.227 

2.227 

2.227 

1.766 

1.766 

3 

1 

3 

444.112 

6.857 

1.930 

1.758 

1.930 

6.618 

2.231 

2.231 

2.231 

1.758 

1.758 

4 

2 

1 

452.937 

6.758 

2.245 

2.019 

2.245 

6.292 

2.116 

2.116 

2.116 

2.019 

2.019 

5 

2 

2 

454.974 

7.177 

1.752 

1.620 

1.752 

7.206 

2.407 

2.407 

2.407 

1.620 

1.620 

6 

2 

3 

443.173 

6.917 

1.782 

1.638 

1.782 

6.903 

2.302 

2.302 

2.302 

1.638 

1.638 

7 

3 

1 

1 263.502 

19.029 

6.261 

5.633 

6.261 

17.377 

5.902 

5.902 

5.902 

5.633 

5.633 

8 

3 

2 

1 282.120 

19.309 

6.354 

5.716 

6.354 

17.633 

5.989 

5.989 

5.989 

5.716 

5.716 

9 

3 

3 

1 641.241 

24.717 

8.133 

7.317 

8.133 

22.572 

7.666 

7.666 

7.666 

7.317 

7.317 

10 

4 

1 

2 035.628 

30.657 

10.088 

9.076 

10.088 

27.996 

9.508 

9.508 

9.508 

9.076 

9.076 

11 

4 

2 

5 853.529 

88.155 

29.008 

26.097 

29.008 

80.504 

27.342 

27.342 

27342 

26.097 

26.097 

12 

4 

3 

756.686 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

13 

5 

4 

460.207 

6.950 

2.282 

2.029 

2.282 

6.346 

2.156 

2.156 

2.156 

2.029 

2.029 

14 

5 

5 

2 035.628 

30.657 

10.088 

9.076 

10.088 

27.996 

9.508 

9.508 

9.508 

9.076 

9.076 

15 

5 

6 

5 853.529 

88.155 

29.008 

26.097 

29.008 

80.504 

27.342 

27.342 

27.342 

26.097 

26.097 

16 

5 

7 

756.686 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

5.180 

17 

0 

8 

771.019 

5.279 

5.279 

5.279 

5.279 

5.279 

5.279 

5.279 

5.279 

5.279 

5.279 

SUMT" 

450.454 

7.090 

1.725 

1.592 

1.725 

7.188 

2.388 

2.388 

2.388 

1.592 

1.592 


‘Sequence of Unconstrained Minimizations Technique. 


TABLE 21. -60-BAR TRUSSED RING: 
LOAD SPECIFICATIONS 


Problem 

Load 

conditions 

Load components, 
kips 



Node 

Px 

p y 


1 

1 

-10.0 

0.0 

Bar60.a 

7 

9.0 

0.0 

Bar60.b 

Bar60.c 

II 

15 

-8.0 

3.0 

Bar60.d 

18 

-8.0 

3.0 


III 

22 

-20.0 

10.0 

Bar60.c 

Lumped 

4 

m x - m y 

, = 200 lb 

Bar60.e 

masses 

12 

m x ~ m y 

, = 100 lb 

Bar60.f 








TABLE 22, -60-BAR TRUSSED RING: 
CONSTRAINT SPECIFICATIONS 


Problem 

Constraint 

type 

Constraint description 

Bar60.a 

Stress 

O/S o {) (i = 1 , 2 ,..., 60) 
Co » 10 ksi 


Minimum area 

Xi = 0.5 in. 2 (/ = 1,2 25) 

Bar60.b 

Displacement 

Constraints along y direction 
1.75 in. at node 4 


Minimum area 

Same as Bar60.a 

Bar60.c 

Frequency 
Minimum area 

/^/o;/o=13Hz 
Same as Bar60.a 


Stress 

Same as Bar 60. a 

Bar60.d 

Displacement 

Constraints along both x and y directions 

1.25 in. at node 10 

1.75 in. at node 4 

2.75 in. at node 1 9 

2.25 in. at node 13 


Minimum area 

Same as Bar60,a 

Bar60.e 

Displacement 

Constraints along both x and y directions 

1.25 in. at node 10 

1.75 in. at node 4 

2.75 in. at node 1 9 

2.25 in. at node 13 


Frequency 

Same as Bar60.c 


Minimum area 

Same as Bar60.a 


Stress 

Same as Bar60.a 

Bar60.f 

Displacement 

Constraints along y direction 
150 in. at node 4 


Frequency 

Same as BarbO.c 


Minimum area 

Same as Bar60.a 


SUMT results. Methods 4, 5, 7, 8, and 9 have large differ- 
ences. For problem Bar60.a (stress constraints only), the fully 
utilized design method, as well as all hybrid methods, per- 
forms well. Performance of the other methods for Bar60,a is, 
in general, inadequate. For problem Bar60.d (stress and dis- 
placement constraints), performance of method 17 is accept- 
able with a 5 percent difference in comparison to the SUMT 
results, method 5 shows less than 10 percent difference, but 
the performance is poor for other methods. For problem 
Bar60.f (all three constraint types), method 5 gives a design 
that is within 5 percent of the SUMT design. The next best 
method for this last problem is method 13, with a 20 percent 
heavier design than SUMT. 


TABLE 23. -60-BAR TRUSSED RING: 
PROBLEMS BAR60.a, BAR60.b, 
BAR60.C, BAR60.d, BAR60.e, 
AND BAR60.f DESIGN 
VARIABLE LINKAGE 


Serial 

number 

Design 

variable 

Members linked 

1 

1 

49,50,51 ,52,53,54, 




2 

2 

1,13 

3 

3 

2,14 

4 

4 

3,15 

5 

5 

4,16 

6 

6 

5,17 

7 

7 

6,18 

8 

8 

7,19 

9 

9 

8,20 

10 


9,21 

11 

11 


12 

12 

11,23 

13 

13 

12,24 

14 

14 

25,37 

15 

15 

26,38 

16 

16 

27,39 

17 

17 


18 

18 

29,41 1 

19 

19 


20 

20 

31,43 

21 

21 

32,44 

22 

22 

33,45 

23 

23 

34,46 

24 

24 

35,47 

25 

25 

36,48 


£. Intermediate Complexity Wing 

The Intermediate Complexity Wing (IC Wing) shown in 
figure 9 is considered next. The finite element model of the 
IC Wing has 88 grid points and a total of 158 elements con- 
sisting of 39 bars, 2 triangular membranes, 62 quadrilateral 
membranes, and 55 shear panels. The wing is made of alumi- 
num and subjected to two load conditions as given in 
table 25. The constraints of the problem are given in table 26. 
The 158 elements are linked to obtain a reduced set of 
57 design variables for optimization as given in table 27. 
Three desTgn cases are considered: (1) IC-Wing.a, for stress 
constraints; (2) IC-Wing.b, for displacement constraints; and 
(3) IC-Wing.c, for stress and displacement constraints. For 
problem IC-Wing.a, the SUMT solution shows 108 active 
stress constraints. For problem IC-Wing.b, it produces three 
active displacement constraints. For IC-Wing.c, no displace- 
ment constraint becomes active, but there are 108 active 
stress constraints. 

The results obtained for the three design cases by 14 of the 
OC-type methods, normalized by SUMT, are presented in 
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TABLE 24.- 60-BAR TRUSSED RING; LINKED DESIGN CONFIGURATION 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Normalized weight 


X 

X 

Bar60.a 

Bar60.b 

Bar60.c 

Bar60.d 

Bar60.e 

Bar60.f 




(stress) 

(displace- 

(frequency) 

(stress 








ment) 


and 











displace- 











ment) 



1 

1 

1 

1.59 

1.05 

1.04 

1.48 

1.03 

1.26 

2 

1 

2 

2.24 

1.00 

1.02 

1.42 

1.01 

2.73 

3 

1 

3 

1.80 

1.00 

1.00 

1.76 

1.00 

3.08 

4 

2 

1 

2.16 

1.05 

(a) 

1.69 

Failed 

1.48 

5 

2 

2 

1.30 

1.03 

1.19 

1.09 

1.14 

1.05 

6 

2 

3 

(a) 

1.01 

(a) 

(a) 

1.08 

1.45 

7 

3 

1 

4.63 

1.00 

1.06 

2.44 

1.76 

3.03 

8 

3 

2 

1.81 



1.02 

2.36 

1.32 

1.39 

9 

3 

3 

(a) 



1.00 

(a) 

1.43 

Failed 

10 

4 

1 

(a) 



1.06 

2.29 

1.08 

4.04 

11 

4 

2 

2.60 



1.02 

3.96 

1.07 

4.09 

12 

4 

3 

(a) 



1.00 

1.63 

1.05 

1.32 

13 

5 

4 

2.11 



2.14 

1.39 

1.02 

1.20 

14 

5 

5 

1.00 



1.06 

1.30 

1.08 

Failed 

15 

5 

6 





1.02 

1.29 

1.07 

Failed 

16 

5 

7 



1 

r 

1.00 

1.29 

1.05 

1.32 

17 

0 

8 



1.33 

1.61 

1.05 

1.61 

1.54 

SUMT* 


1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

Weight, lb 

286.862 

294.066 

391.902 

308.730 

391.820 

420.243 


“Method failed for this problem. 

b Sequence of Unconstrained Minimizations Technique. 



Figure 9. — Intermediate complexity wing. (Representative elements are circled, 
nodes are not.) 
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TABLE 25. -INTERMEDIATE COMPLEXITY WING: LOAD SPECIFICATIONS 
(a) Load condition I 


Load components, 
lb 

Node 


p y 

p z 

53 

0 

0 

404.0 

54 





404.0 

55 





420.0 

56 

1 


1 


420.0 

57 

-4640.0 

1900.0 

1120.0 

58 

4640.0 

-1900.0 

1120.0 

59 

2290.0 

-937.0 

883.0 

60 

-2290.0 

937.0 

883.0 

61 

0 

0 

413.0 

62 





413.0 

63 





391.0 

64 





391.0 

65 





368.0 

66 


f 


r 

368.0 

67 

-3030.0 

1240.0 

804.0 

68 

3030.0 

-1240.0 

804.0 

69 

3070.0 

-520.0 

1040.0 

70 

-3070.0 

520.0 

1040.0 

71 

0 

0 

433.0 

72 





433.0 

73 





370.0 

74 





370.0 

75 





304.0 

76 





304.0 

77 

-1370.0 

262.0 

446.0 

78 

1370.0 

-262.0 

446.0 


Load components, 
lb 

Node 

p . 

p y 


27 

-4070.0 

1660.0 

902.0 

28 

4070.0 

-1660.0 

902.0 

29 

1740.0 

-713.0 

646.0 

30 

-1740.0 

713.0 

646.0 

31 

0 

0 

340.0 

32 





340.0 

33 





352.0 

34 





352.0 

35 





365.0 

36 





365.0 

37 

-4250.0 

1740.0 

974.0 

38 

4250.0 

-1740.0 

974.0 

39 

1820.0 

-743.0 

694.0 

40 

-1820.0 

743.0 

694.0 

41 

0 

0 

365.0 

42 





365.0 

43 





378.0 

44 





378.0 

45 





392.0 

46 





392.0 

47 ! 

-4440.0 

1820.0 

1050.0 

48 ; 

4440.0 

-1820.0 

1050.0 

49 

1890.0 

-773.0 

742.0 

50 

-1890.0 

773.0 

742.0 

51 



0 

390.0 

52 


3 

0 

390.0 


(b) Load condition II 


Load components, 

lb 

Node 

p . 

p y 

P z 

1 

205.0 

7380.0 

926.0 

2 

-205.0 

7380.0 

926.0 

3 

0 

0 

29.0 

4 

0 

0 

29.0 

5 

-2800.0 

-6960.0 

1130.0 

6 

2800.0 

6960.0 

1130.0 

7 

0 

0 

90.9 

8 

0 

0 

90.9 

9 

-9870.0 

-9780.0 

1130.0 

10 

9870.0 

9780.0 

1130.0 

11 

0 

0 

178.0 

12 





178.0 

13 





214.0 

14 





214.0 

15 





253.0 

16 





253.0 

17 

-5680.0 

2320.0 

1020.0 

18 

5680.0 

2320.0 

1020.0 

19 

2310.0 

-946.0 

723.0 

20 

-2310.0 

946.0 

723.0 

21 

0 

0 

314.0 

22 





314.0 

23 





326.0 

24 





326.0 

25 





338.0 

26 





338.0 


Load components, 
lb 

Node 

p . 


P z 

53 

0 

0 

334.0 

54 





334.0 

55 





264.0 

56 



i 


264.0 

57 

-1380.0 

565.0 

386.0 

58 

1380.0 

-565.0 

386.0 

59 

5300.0 

-2170.0 

1820.0 

60 

-5300.0 

2170.0 

1820.0 

61 

0 

0 

458.0 

62 





458.0 

63 





326.0 

64 





326.0 

65 





233.0 

66 

i 




233.0 

67 

-922.0 

377.0 

287.0 

68 

922.0 

-377.0 

287.0 

69 

7160.0 

-1210.0 

2180.0 

70 

-7160.0 

-1210.0 

2180.0 

71 


) 

0 

484.0 

72 





484.0 

73 





310.0 

74 





310.0 

75 





194.0 

76 





194.0 

77 

-451.0 

86.0 

175.0 

78 

451.0 

-86.0 

175.0 


Load components, 
lb 

Node 

r. 


P z 

27 

1210.0 

496.0 

311.0 

28 

1210.0 

—496.0 

311.0 

29 

3990.0 

-1630.0 

1310.0 

30 

-3990.0 

1630.0 

1310.0 

31 

0 

0 

375.0 

32 





375.0 

33 





291.0 

34 





291.0 

35 





230.0 

36 




r 

230.0 

37 

-1270.0 

518.0 

336.0 

38 

1270.0 

-518.0 

336.0 

39 

4160.0 

-1700.0 

1410.0 

40 

-4160.0 

1700.0 

1410.0 

41 

0 

0 

402.0 

42 





402.0 

43 





313.0 

44 





313.0 

45 





247.0 

46 




f 

247.0 

47 

-1320.0 

541.0 

361.0 

48 

1320.0 

-541.0 

361.0 

49 

4330.0 

-1770.0 

1500.0 

50 

-4330.0 

1770,0 

1500.0 

51 

< 

) 

0 

430.0 

52 

!L 



D 

430.0 


Load components, 
lb 

Node 

p z 

p y 

p z 

1 

351.0 

-12 600.0 

1530.0 

2 

-351.0 

12 600.0 

1530.0 

3 

0 

0 

29.5 

4 

0 

0 

29.5 

5 

-2420.0 

-6020.0 

979.0 

6 

2420.0 

6020.0 

979.0 

7 

0 

0 

55.9 

8 

0 

0 

55.9 

9 

-4020.0 

-3980.0 

474.0 

10 

4020.0 

3980.0 

474.0 

11 

0 

0 

194.0 

12 





194.0 

13 





175.0 

14 





175.0 

15 





157.0 

16 


r 


r 

157.0 

17 

-1600.0 

653.0 

325.0 

18 

1600.0 

-653.0 

325.0 

19 

5510.0 

-2250.0 

1550.0 

20 

-5510.0 

2250.0 

1550.0 

21 

0 

0 

347.0 

22 





347.0 

23 





270.0 

24 





270.0 

25 





213.0 

26 





213.0 


24 
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TABLE 26. - INTERMEDIATE COMPLEXITY (IC) WING: 
CONSTRAINT SPECIFICATIONS 


Problem 

Constrain! 

type 

Constraint description 

IC-Wing.a 

Stress 

Minimum area 

0,^00 (/= 1,2,...,158) 

o 0 = 10.5 ksi 

X; = 0.001 in. 2 (f = 1,2,. .,57) 

IC-Wing.b 

Displacement 

Constraints along traverse direction 
10 in. at node 1 
10 in. at node 10 


Minimum area 

Same as IC-Wing.a 


Stress 

Same as IC-Wing.a 

IC-Wing.c 

Displacement 

Same as IC-Wing.b 


Minimum area 

Same as IC-Wing.a 


TABLE 27. - INTERMEDIATE COMPLEXITY (IC) WING: PROBLEMS IC-WING.a, IC-WING.b, 
AND IC-WING.c DESIGN VARIABLE LINKAGE 


Serial 

number 

Design 

variable 

Members linked 

1 

1 

1,2 

2 

2 

3,4 

3 

3 

5,6 

4 

4 

7,8 

5 

5 

9,10 

6 

6 

11,12 

7 

7 

13,14 

8 

8 

15,16 

9 

9 

17,18 

10 

10 

19,20 

11 

11 

21,22 

12 

12 

23,24 

13 

13 

25,26 

14 

14 

27,28 

15 

15 

29,30 

16 

16 

31,32 

17 

17 

33,34 

18 

18 

35,36 

19 

19 

37,38 

20 

20 

39,40 

21 

21 

41,42 

22 

22 

43,44 

23 

23 

45,46 

24 

24 

47,48 

25 

25 

49,50 

26 

26 

51,52 

27 

27 

53,54 

28 

28 

55,56 

29 

29 

57,58 

30 

30 

59,60 

31 

31 

61,62 


Serial 

number 

Design 

variable 

Members linked 

32 

32 

63,64 

33 

33 

65,66,67,68,69,70,71,72,73,74,75, 
76,77,78,79,80,81,82,83,84,85,86, 
87,88,89,90,91 ,92,93,94,95,96 

34 

34 

97 

35 

35 

98 

36 

36 

99 

37 

37 

100 

38 

38 

101 

39 

39 

102 

40 

40 

103 

41 

41 

104 

42 

42 

105 

43 

43 

106 

44 

44 

107 

45 

45 

108 

46 

46 

109 

47 

47 

110 

48 

48 

111 

49 

49 

112 

50 

50 

113 

51 

51 

114 

52 

52 

115 

53 

53 

116 

54 

54 

117 

55 

55 

118 

56 

56 

119 

57 

57 

120,121,122,123,124,125,126,127,128,129, 

130,131,132,133,134,135,136,137,138,139, 

140,141,142,143,144,145,146,147,148,149, 

150,151,152,153,154,155,156,157,158 





TABLE 28.- INTERMEDIATE COMPLEXITY (IC) WING: 
LINKED DESIGN CONFIGURATION 


Methods 

Optimality 

criteria 

(O Q 

formulas 

Normalized weight 


X 

X 

IC-Wing.a 

IC-Wing.b 

IC-Wing.c 




(stress) 

(displace- 

(stress and 






ment) 

displacement) 

1 

1 

1 

(a) 

(a) 

( a ) 

2 

1 

2 

(a) 

(a) 

(a) 

3 

1 

3 

(a) 

(a) 

(a) 

4 

2 

1 

(b) 

.99 

(b) 

5 

2 

2 

3.22 

1.01 



6 

2 

3 

(b) 

1.11 



7 

3 

1 

(b) 

1.05 



8 

3 

2 

1.53 

1.05 



9 

3 

3 

3.93 

1.16 



10 

4 

1 

(b) 

1.06 



11 

4 

2 

(b) 

1.06 



12 

4 

3 

0>) 

1.69 



13 

5 

4 

2.93 

1.06 



14 

5 

5 

1.00 

1.06 



15 

5 

6 



1.06 



16 

5 

7 



1.69 



17 

0 

8 


f 

1.63 

LOO 

1 SUMT 0 


1.00 

LOO 

LOO 

Weight, lb 

388.390 

84.114 

391.902 


'Methods I to 3 arc not applicable. 

^Method failed for this problem. 

Sequence of Unconstrained Minimizations Technique. 


table 28. Note that the linear form of the Lagrange multiplier 
update formula is not applicable for element types with mul- 
tiaxial stress fields, such as membranes. For stress constraints 
only (IC-Wing.a), only the fully utilized design method and 
the hybrid methods perform satisfactorily. For problem 
IC-Wingx, where displacement constraints were added but 
did not become active in the SUMT solution, even the hybrid 
methods fail. For all methods in which matrix inversion of a 
constraint sensitivity matrix is required, matrix singularities 
are detected during inversion. Floating point divide excep- 
tions are encountered in most of the other methods for this 
problem. Only the fully utilized design method is successful 
in this stress and displacement problem. For displacement 
constraints only (IC-Wing.b), nine of the methods (4, 5, 7, 8, 
10, 11, 13, 14, and 15) perform well, two (6 and 9) are mar- 
ginal, and two (12 and 16) are poor. 

F. Forward-Swept Wing 

The finite element model of the forward-swept wing 
(FSW), made of aluminum and shown in figure 10, has 
30 grid points and a total of 135 truss elements. The loads are 
given in table 29. The constraints of the problem are given in 
table 30. Two design cases are considered. The first one treats 
all 135 element areas as independent design variables, 
whereas in the second case the bar areas are linked to obtain 
a reduced set of 61 linked design variables, as shown in 
table 31. For each case, three situations are considered: sets 
FSW.a, FSW.b, and FSW.c include stress, displacement, and 
frequency constraints, respectively, for the unlinked problem; 
whereas sets FSW.d, FSW.e, and FSW.f include stress, dis- 
placement, and frequency constraints, respectively, for the 
linked problem. The number of active constraints found by 



Figure 10. — Forward-swept wing. (Representative elements are circled, nodes are not.) 
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TABLE 29. -OPTIMUM DESIGN OF FORWARD SWEPT WING: 
LOAD SPECIFICATIONS 


Problem 

Load 

conditions 

Load components, kips 

Node 

p. 

Py 

Py 

FSW.a 

FSW.b 

FSW.e 

FSW.d 

FSW.e 

FSW.f 

■ 

9 

10 

29 

30 

0.0 

0.0 

40.0 

40.0 

20.0 
20.0 

FSW.e 

FSW.f 

Lumped 

masses 

(1,2,. ..,8,10) 

9 

(1 1,1 2,..., 20) 

(21, 22, ...,25,27, 28,. ..,30) 
26 

m x ~ rriy = m z = 
n\ x — n\y — m z 

ttt x — tfiy — tti z 
fft x = Tf\y — TTl z 
m x ~ fTiy — m z 

= 151b 
= 100 lb 
= 30 lb 
= 18 lb 
= 300 lb 


TABLE 30. -FORWARD-SWEPT WING (FSW): 
CONSTRAINT SPECIFICATIONS 


Problem 

Constraint 

Constraint description 


type 



Stress 

a ( £a 0 (/ = 1,2, ...,135) 

FSW.a 


a 0 = 10 ksi 

FSW.d 

Minimum area 

Xi = 0.25 in. 2 (i = 1 35) for FSW.a 
Xj = 0.25 in. 2 (i = 1,2 61) for FSW.d 

FSW.b 

FSW.e 

Displacement 

Constraints along traverse direction 
10.0 in. at node 10 
10.0 in, at node 30 


Minimum area 

Same as FSW.a or FSW.d 

FSW.e 

FSW.f 

Frequency 

/^/ u; / 0 =10Hz 


Minimum area 

Same as FSW.a or FSW.d 


TABLE 31. -FORWARD-SWEPT WING (FSW): PROBLEMS FSW.d, FSW.e, 
AND FSW.f DESIGN VARIABLE LINKAGE 


Serial 

number 

Design 

variable 

Members linked 

1 

1 

121,122,123,124,125, 



126,127,128,129,130, 



131,132,133,134,135 

2 

2 

1,2 

3 

3 

3,4 

4 

4 

5,6 

5 

5 

7,8 

6 

6 

9,10 

7 

7 

11,12 

8 

8 

13,14 

9 

9 

15,16 

10 

10 

17,18 

11 

11 

19,20 

12 

12 

21,22 

13 

13 

23,24 

14 

14 

25,26 

15 

15 

27,28 

16 

16 

29,30 

17 

17 

31,32 

18 

18 

33,34 

19 

19 

35,36 

20 

20 

37,38 

21 

21 

39,40 

22 

22 

41,42 

23 

23 

43,44 

24 

24 

45,46 

25 

25 

47,48 

26 

26 

49,50 

27 

27 

51,52 

28 

28 

53,54 

29 

29 

55,56 

30 

30 

57,58 


Serial 

number 

Design 

variable 

Members linked 

31 

31 

59,60 

32 

32 

61,62 

33 

33 

63,64 

34 

34 

65,66 

35 

35 

67,68 

36 

36 

69,70 

37 

37 

71,72 

38 

38 

73,74 

39 

39 

75,76 

40 

40 

77,78 

41 

41 

79,80 

42 

42 

81,82 

43 

43 

83,84 

44 

44 

85,86, 

45 

45 

87,88 

46 

46 

89,90 

47 

47 

91,92 

48 

48 

93,94 

49 

49 

95,96 

50 

50 

97,98 

51 

51 

99,100 

52 

52 

101,102 

53 

53 

103,104 

54 

54 

105,106 

55 

55 

107,108 

56 

56 

109,110 

57 

57 

111,112 

58 

58 

113,114 

59 

59 

115,116 

60 

60 

117,118 

61 

61 

119,120 


27 




/ 






SUMT for the six problems (FSW.a to FSW.f) are 75 stress, 
2 displacement, 1 frequency, 68 stress, 2 displacement, and 
1 frequency constraint, respectively. 

Results for the six problems solved by the 17 methods are 
given in tables 32 and 33. Many of the first 16 methods per- 
form adequately for problems FSW.b, FSW.c, FSW.e, and 
FSW.f (displacement or frequency constraints). Most of the 
methods that include a reciprocal form for the design variable 
update formula (3, 9, 12, and 16) perform very marginally, at 
best, for problems with displacement constraints (FSW.b and 
FSW.e). The only other method that includes a reciprocal 
form (method 6) has some difficulty with frequency con- 
straints (FSW.c and FSW.f). Method 13 does not perform 
well on the frequency-constrained problems either. For prob- 
lems FSW.a and FSW.d (stress constraints), only the hybrid 
methods and the fully utilized design methods perform well, 
with differences within 1 percent of the SUMT results. 

G. Three-Bar Truss 

Several additional features are illustrated with a three-bar 
truss, as shown in figure 1. The first two problems are con- 


TABLE 32. -FORWARD SWEPT WING (FSW) 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Normalized weight 

X 

X 

FSW.a 

(stress) 

FSW.b 

(displace- 

ment) 

FSW.c 

(frequency) 

1 

1 

1 

5.64 

1.19 

1.09 

2 

1 

2 

5.39 

1.15 

1.08 

3 

1 

3 

(a) 

2.37 

1.07 

4 

2 

1 

9.39 

1.09 

1.00 

5 

2 

2 

1.31 

1.05 

1.00 

6 

2 

3 

1.67 

1.03 

1.18 

7 

3 

1 

19.47 

1.09 

1.01 

8 

3 

2 

1.53 

1.09 



9 

3 

3 

1.98 

1.40 



10 

4 

1 

8.36 

1.07 



11 

4 

2 

1.57 

1.07 



12 

4 

3 

(a) 

1.33 


r 

13 

5 

4 

1.16 

1.08 

3,01 

14 

5 

5 

1.01 

1.07 

1,01 

15 

5 

6 



1.07 

1.01 

16 

5 

7 



1.33 

1.01 

17 

0 

8 



2.80 

2.76 

SUMT 11 


1.00 

1.00 

1.00 

Weight, lb 

2 793.179 

672.574 

230.045 


“Method failed for this problem. 

b Sequence of Unconstrained Minimizations Technique. 


TABLE 33. - FORWARD SWEPT WING (FSW) 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Normalized weight 

X 

X 

FSW.d 

(stress) 

FSW.e 

(displace- 

ment) 

FSW.f 

(frequency) 

1 

1 

1 

5.23 

1.17 

1.10 

2 

1 

2 

3.72 

1.13 

1.09 

3 

1 

3 

(a) 

2,34 

1.07 

4 

2 

1 

5.64 

1.08 

1.00 

5 

2 

2 

1.24 

1.05 

1.00 

6 

2 

3 

1.32 

1.03 

1.19 

7 

3 

1 

24.15 

1.10 

1.02 

8 

3 

2 

1.58 

1.10 



9 

3 

3 

2.26 

1.39 



10 

4 

1 

19,76 

1.07 



11 

4 

2 

1.49 

1.07 



12 

4 

3 

(a) 

1.40 


r 

13 

5 

4 

1.60 

1.07 

2.99 

14 

5 

5 

1.00 

1.07 

1.02 

15 

5 

6 



1.07 

1.02 

16 

5 

7 



1.40 

1.02 

17 

0 

8 

1 


2.73 

2.74 

SUMT* 


1.00 

1.00 

LOO 

Weight, lb 

3 009.874 

690.338 

231.803 


“Method failed for this problem. 

^Sequence of Unconstrained Minimizations Technique. 


TABLE 34. -THREE-BAR TRUSS: 
LOAD SPECIFICATIONS 


Problem 

Load 

conditions 

Load components, 
kips 



Node 

/’x 

r, 


I 

1 

50.0 

100.0 

Bar3.a 

11 

1 

-50.0 

-100.0 


II! 

1 

50.0 

0.0 

Bar3.b 

i 

1 

50.0 

100.0 


, 

1 

70.0 

0.0 

Bar3.c 

ii 

1 

-35.0 

-95.0 

Bar3.a 

Bar3.b 

Lumped 

masses 

1 

fft X — TTly 

= 262.5 lb 
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TABLE 35.— THREE-BAR TRUSS; 
CONSTRAINT SPECIFICATIONS 


Problem 

Constraint 

type 

Constraint description 


Stress 

a, <ct 0 (r = 1,2,3) 
CTo = 20 ksi 

Bar3.a 

Bar3.b 

Displacement 

Constraints at node 1 of magnitude: 
0.20 in. along the x direction 
0.05 in. along the y direction 


Frequency 

/>/ 0 ;/ 0 = 114.3 Hz 


Minimum area 

x, = 0.001 in. 2 (r = 1,2,3) 

Bar3.c 

Stress 

Gi<c 0 (i =1,2,3) 
— 15 ksi 


Minimum area 

Same as Bar3.a 


sidered with three load conditions (for Bar3.a) and a single 
load condition (for Bar3.b), as given in table 34. Note that the 
load condition in problem Bar3.b is identical to the first load 
condition in problem Bar3.a. The stress, displacement, and 
frequency constraints are specified in table 35. SUMT pro- 
vides an optimal solution with one active frequency, two 
active displacement, and three active stress constraints for 
problem Bar3.a, as well as one active frequency, one active 
displacement, and one active stress constraint for problem 
Bar3.b. The results of these two problems, shown in tables 36 
and 37, are discussed in section V. 

The third problem in this set (Bar3.c) considers six cases 
with stress constraints under two load conditions (see 
tables 34 and 35). In the first case, the densities of the three 
bars are taken to be the same (0.1 lb/in. 3 ). In each of the sub- 
sequent cases, the densities are allowed to differ. The six 
cases within this problem are each solved with several 
optimizers, including OC (refs. 3 and 4), SUMT (ref. 21), the 
International Mathematical and Statistical Library’s 


TABLE 36 — OPTIMUM DESIGN OF THREE-BAR TRUSS: PROBLEM BAR3.a 
WITH THREE LOAD CONDITIONS 

[Stress, displacement, and frequency constraints. Case 1 : six active constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design 
variables, 
in. 2 

Active Constraints 

X 

X 

Xi 

X2 

X3 

Stress 

Displacement 

Frequency 

1 

1 

1 

118.475- 

1.993 

3.863 

3.652 

(b) 

S 1 1 *6 1 

(b) 

2 

1 

2 

118.287 

1.985 

3.866 

3.646 

(b) 

g]l»gl3 



3 

l 

3 

116.807 

1.917 

3.882 

3.598 

g3>g6 

gn>gi3 



4 

2 

1 

166.293 

4.247 

2.852 

5.495 

(b) 

gll*gl3 


r 

5 

2 

2 

109.423 

1.008 

4.726 

3.388 


<b) 

gl6 

6 

2 

3 

128.181 

2.423 

3.700 

4.025 

(b) 

gn»gi3 

(b) 

7 

3 

1 

15 460.000 

183.667 

535.264 

531.058 



(b) 

(b) 

8 

3 

2 

15 948.900 

191.988 

566.978 

534.860 



(b) 

(b) 

9 

3 

3 

147 900.000 

(a) 

271.077 

271.077 



(b) 

gl6 

10 

4 

1 

123.696 

1.080 

3.166 

5.428 


r 

gll»gl3 

(b) 

11 

4 

2 

105.878 

1.224 

3.718 

3.634 

g3’g6>g? 

gll’gl3 

gl6 

12 

4 

3 

107.537 

1.251 

3.671 

3.757 

(b) 

gl!>gl3 

(b) 

13 

5 

4 

117.430 

1.445 

3.432 

4.432 



gtl»g]3 



14 

5 

5 

126.122 

1.141 

3.128 

5.565 



gll’g]3 



15 

5 

6 

205.515 

6.429 

1.939 

6.732 



glbgl3 



16 

5 

7 

14 390.000 

.151 

360.764 

762.766 



(b) 



17 

0 

8 

122.221 

1.571 

3.333 

4.714 



glbg|3 



SUMT 

100.078 

1.091 

3.848 

3.264 

g3>g6»g7 

gtbgl3 

g 16 


a Area greater than 1000 in 2 . 
b No active constraints of this type. 

Sequence of Unconstrained Minimizations Technique. 
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TABLE 37.— OPTIMUM DESIGN OF THREE-BAR TRUSS: PROBLEM BAR3.b WITH ONE 

LOAD CONDITION 

[Stress, displacement, and frequency constraints. Case 2: three active constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Optimum design 
variables, 
in. 2 

Active constraints 

X 

X 

Xi 

X2 

X3 

Stress 

Displacement 

Frequency 

1 

1 

1 

115.125 

1,840 

3.899 

3.543 

g3 

g5 

(b) 

2 

1 

2 

117.322 

1.940 

3.876 

3.615 

(b) 

g5 

(b) 

3 

1 

3 

117.598 

1.953 

3.874 

3.624 

(b) 

g5 

(b) 

4 

2 

1 

100.085 

1.091 

3,851 

3.262 

g3 

g5 

g6 

5 

2 

2 

112.170 

1.470 

3.661 

3.873 

(b) 

g5 

(b) 

6 

2 

3 

150.214 

.001 

6.070 

6.373 



(b) 

g6 

7 

3 

1 

13 000.000 

242.722 

274.612 

480.441 





(b) 

8 

3 

2 

14 300.000 

142.087 

583.845 

455.524 







9 

3 

3 

12 400.000 

386.746 

0.000 

494.501 







10 

4 

1 

53 800.000 

.001 

2 232.740 

2 232.470 







11 

4 

2 

46 200.000 

1 206.536 

1 206.536 

1 206.536 






f 

12 

4 

3 

147 000.000 

270.881 

270.881 

(a) 





g6 

13 

5 

4 

135.991 

.257 

4.613 

6.097 





g6 

14 

5 

5 

180.738 

2.931 

3.037 

7.701 





(b) 

15 

5 

6 

4 444.091 

72.560 

155.956 

131.407 





(b) 

16 

5 

7 

611.774 

36.218 

3.678 

4.440 



g5 

(b) 

17 

0 

8 

150.728 

.002 

5.023 

7.104 


r 

(b) 

g6 

SUMT 

100.074 

1.089 

3.848 

3.266 

g3 

g5 

g6 


a Area greater than 1000 in 2 . 
b No active constraints of this type. 

Sequence of Unconstrained Minimizations Technique. 


nonlinear optimization routine based on sequential quadratic 
programming (IMSL) (ref. 9), another sequential quadratic 
programming code (SQP) (ref. 7), a sequential linear 
programming code (SLP) (ref. 8), and feasible directions 
(FD) (ref. 8). These results are summarized in table 38 and 
examined in section V. 

V. Discussion 

In this section the overall performance of the methods are 
examined on the basis of the numerical examples solved, and 
some related topics are discussed,, 

A. Performance of Optimality Criteria Methods 

The derivable OC method was developed first for displace- 
ment constraints and was later extended for stiffness con- 
straints that can be stated in terms of work quantities. The 
external work of a virtual load system along the constrained 
displacement is an example. Convergence, as in the case of the 


heuristic fully stressed design, is a function only of the sensitiv- 
ity of the internal forces to changes in the number of truss 
members. The method exhibits excellent performance if the 
behavior approximates that of statically determinate structures. 
Convergence is not a function of the number of design vari- 
ables, Buckling and vibration frequency constraints also can be i 
stated in energy terms and brought under the OC formulation, 1 
but because of the nature of the eigenvalue problem, OC § 

convergence can be dependent on the specific nature of the | 

structure. 

The elegancy of the criteria for displacement constraints is 
illustrated by considering a set of four large, three-dimensional 
trusses with 148 to 1027 design variables (figs. 1 1(a) to (d)). 

The trusses shown in figure 1 1 (a) have four concentrated loads 
of different magnitudes (20 and 40 kip) which simultaneously I 
induces twisting and flexing of the cantilevered trussed slabs. 

Two displacement constraints are imposed: the displacement 
must not exceed (1) 10 in. at the higher load point and (2) 20 in. 
at the lower load point (fig. 11(a)), and at optimum both 
displacement constraints are active. The problem is solved 
using the OC method. Design weight versus iteration history 
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TABJ.E 38. — OPTIMUM DESIGNS OFTHREE-BAR TRUSS: PROBLEM BAR3.C 
WITH TWO LOAD CONDITIONS 


[Number of active constraints exceed number of design variables.] 


Method* 

Cost 

coefficients 

Member areas 

Active 

constraints 

Optimum 

weight, 

lb 

Remarks 

SUMT 

0.1 

0.1 0.1 

3.291 

3.986 

3.323 

81 * 83 * 85*86 

L33X10 2 




3 

6 8 

3.299 

3.998 

3.299 



7.53X10* 




6 

12 18 

3.299 

3.997 

3.298 



1.59X10 4 




16 

13 25 

3.298 

3.998 

3.299 


► 

2.43 X10 4 




1 

200 300 

67.068 

9.111 

0.001 

B5>86 

1.92X10 5 

(c) 


1 

5 900 

56.629 

-1.339 

-1.937 

<*>) 

—2.39X10* 

(d) 

OC 14 

.1 

.1 .1 

3.299 

3.999 

3.300 

81 * 83 * 85*86 

1.33X10 2 




3 

6 8 









7.53X10* 




6 

12 18 









1.60X10 4 




16 

13 25 









2.43 X10 4 




1 

200 300 









2.20 X10 5 




1 

5 900 


f 


f 


f 


r 

4.22 X 10 s 



007 

.1 

.1 .1 

3.299 

3.999 

3.300 

81 - 83 * 85*86 

1.33X10 2 




3 

6 8 









7.53X10* 




6 

12 18 









1.60X1 0 4 




16 

13 25 









2.43 X10 4 




1 

200 300 









2.20X1 0 5 




1 

. 5 900 


f 


f 


f 


f 

4.22X10 5 



OC 2 

.1 .1 .1 

5.469 

.007 

12.655 

8s 

2.56 X JO 2 

(c) 


3 

6 8 

3.234 

3.919 

3.446 

81 * 83 * 85*86 

7.62X10* 




6 

12 18 

3.242 

3.936 

3.426 

81 * 83 * 85*86 

1.62X10 4 




16 

13 25 

3.251 

3.958 

3.407 

81 * 83 * 85*86 

2.45 X10 4 




1 

200 300 

4.839 

3.776 

3.574 

85*86 

2.28 X 10 s 




1 

5 900 

3.340 

4.031 

3.288 

81 * 83 * 85*86 

4.21 X 10 s 



IMSL 

.1 

.1 .1 

1.000 

1.000 

1.000 

(b) 

3.83X10* 

<e) 


3 

6 8 

3.299 

4.000 

3.299 

81 * 83 * 85*86 

7.53X10* 




6 

12 18 









1.60X10 4 




16 

13 25 









2.43 X10 4 




1 

200 300 


r 


f 


r 

\ 

r 

2.20 X10 5 




1 

5 900 

13.352 

11.509 

.001 

86 

7.77X10* 

(c) 

SOP 

.1 

.1 .1 

2.335 

2.503 

2.505 

0>) 

9.35 X10 1 




3 

6 8 

2.334 







5.32X10* 




6 

12 18 

2.334 







1.14X10 4 




16 

13 25 

2.335 







1.74X10 4 




1 

200 300 

2.335 




r 


r 

1.57X10 5 

(e) 


1 

5 900 

8.144 

8.729 

.001 

85 

5.64X10* 

(d) 

SLP 

.1 

.1 .1 

3.313 

3.971 

3.323 

81 * 83 * 85*86 

1.33X10 2 




3 

6 8 

3.309 

3.963 

3.334 



7.55XIO* 




6 

12 18 

3.309 

3.962 

3.335 



1.60X10 4 




16 

13 25 

3.308 

3.961 

3.336 



2.44X10 4 

1 



1 

200 300 

3.300 

3.967 

3.328 


F 

2.21X10 5 




1 

5 900 

12.630 

11.726 

.001 

83*86 

7.78X10* 

(d) 

FD 

.1 

.1 .1 

3.199 

2.556 

5,102 

81*85 

1.42X10 2 

(c) 


3 

6 8 

3.501 

2.402 

4.682 

(b) 

8.22X10* 

(d) 


6 

12 18 

3.635 

2353 

4.563 

86 

1.75X10 4 

(c) 


16 

13 25 

3.657 

2316 

4.467 

86 

2.71 X10 4 

(c) 


1 

200 300 

3.684 

2.336 

4.512 

86 

2.39 X 10 s 

(c) 


1 

5 900 

3.167 

2.105 

3.822 

81*83 

4.88X10 5 

(c) 


■Methods defined in appendix B. 
b No active constraints, 
incomplete. 

d Mcthod failed for this problem. 
c FInaI design is infeasible. 



Weight, lb 


(c) Grid, 6x6; members, 529; optimum weight, 14 310 lb. 


(d) Grid, 6x12; members, 1 027; optimum weight, 17 455 lb. 


Figure 1 1 . — Three-dimensional truss problems. 


30x1 0 3 



Figure 1 2. — Convergence curves for three-dimensional truss 
problems. 


(fig. 12) indicates smooth monotonic convergence, and the 
number of iterations required for convergence also is indepen- 
dent of the number of design variables for all four examples. The 
OC method performed very well for all four design problems. 

In this study, the optimality criteria method achieves satis- 
factory solutions, closely matching the optimal design 
obtained with the use of the mathematical programming tech- 
nique SUMT, for problems in which only displacement and/ 
or frequency constraints arise. For each of these problems, at 
least one of the OC methods produces a result close to the 
SUMT answer (usually within 1 percent, always within 
3 percent). The most consistent method for these problems 
is OC method 5 (with the use of the exponential form of the 
Lagrange multiplier update formula and the linearized form 
of the design variable update formula). This OC method 
shows a 19 percent heavier design in the worst case (Bar60.c; 
table 24). 

For problems that contain only a few design variables (up 
to 10) for which stresses are constrained, for which displace- 
ments and/or frequencies may also be constrained, the 
optimality criteria method also appears to be satisfactory. In 
the worst case (Bari 0. a; table 17), the best solution gives a 
7 percent heavier design than SUMT. 


For larger problems that include stress constraints, the 
optimality criteria method may not be as effective. The best 
solution for the best case (Bar60.f; table 24) is 5 percent 
over-designed. For the worst case (IC-Wing.c), the solution 
diverged for all of the OC methods (see table 28). The cause 
of the divergence may, however, be due to linear functional 
dependence among active constraints (see ref. 16). For a 
problem with a large number of active constraints that 
include stress limitations, the optimality criteria method 
appears to follow a subset (or subsets) of these active con- 
straints, and possibly for this reason yields a heavier design. 
For example, the optimum solution obtained by SUMT for 
problem Bar60.f shows a total of 19 active constraints, con- 
sisting of 1 frequency, 1 displacement, and 17 stress con- 
straints. The designs of the ring given by the OC methods 
yield over-designs with fewer active constraints (see 
table 39). For example, method 1, which over-designs by 
26 percent, has only six active stress constraints and neither 
displacement nor frequency constraints become active. The 
best solution (method 5) yields a 5 percent over-design with 
nine active stress constraints, as well as active displacement 
and frequency constraints. 

This observation is seen in all the larger problems in 
which stress constraints play a role (see table 40). No OC 


TABLE 39.— ACTIVE CONSTRAINTS BY TYPE AND NORMALIZED WEIGHTS, 
60 BAR TRUSSED RING: PROBLEM BAR60.f 
[Stress, displacement, and frequency constraints.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Normalized 

weight 

Index of active constraints 
by type 

X 

X 

Stress 

Displacement 

Frequency 

1 

1 

1 

1.26 

1,15,31,36,43,48 

(c) 

(c) 

2 

1 

2 

2.73 

36,48 



(c) 

3 

l 

3 

3.08 

36,48 



(C) 

4 

2 

1 

1.48 

U3 



184 

5 

2 

2 

1.05 

1,13,26,30,31,36,42,43,48 

181 



6 

2 

3 

1.45 

(c) 

(c) 



7 

3 

1 

3.03 

6,18 





8 

3 

2 

1.39 

80 





9 

3 

3 

(a) 

(c) 





10 

4 

1 

4.04 

(c) 




f 

11 

4 

2 

4.09 

(c) 



(C) 

12 

4 

3 

1.32 

25,30,36,37,42,48,80 



(c) 

13 

5 

j 4 

1.20 

25,30,36,37,42,48,80 


f 

184 

SIJMT 6 


1.00 

1,7,12,13,25,26,30,31,36, 

181 

184 






37,42,43,48,49,55,80,114 


1 


1 


a Method failed for this problem. 

^Sequence of Unconstrained Minimizations Technique. 
c No active constraints of this type. 
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TABLE 40. - NUMBER OF ACTIVE CONSTRAINTS AND NORMALIZED WEIGHT 
[Active constraint types are stress a, displacement X, and frequency /.] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Bar60.a 

Bar60.d 

Bar60.f 

IC-Wing.a 

FSW.a 

FSW. 

i 

d 


X 

X 

Weight 

cr 

Weight 

a 

X 

Weight 

a 

X 

/ 

Weight 

cr 

Weight 

cr 

Weight 

o 

1 

1 

1 

1.59 

*8 

1.48 

8 

0 

1.26 

6 

0 

0 

(b) 

(b) 

5.64 

2 

5.23 

2 

2 

1 

2 

2.24 

2 

1.42 

9 

0 

2.73 

2 

0 

0 

(b) 

(b) 

5.39 

1 

3.72 

4 

3 

1 

3 

1.80 

6 

1.76 

4 

0 

3.08 

2 

0 

0 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

4 

2 

1 

2.16 

2 

1.69 

2 

0 

1.48 

2 

0 

1 

<b) 

(b) 

9.39 

2 

5.64 

5 

5 

2 

2 

1.30 

15 

1.09 

17 

1 

1.05 

9 

1 

1 

3.22 

1 

1.31 

14 

1.24 

21 

6 

2 

3 

(b) 

(b) 

(b) 

(b) 

(b) 

1.45 

0 

0 

1 

(b) 

(b) 

1.67 

26 

1.32 

19 

7 

3 

1 

4.63 

5 

2,44 

2 

0 

3.03 

2 

0 

1 

(b) 

(b) 

19.47 

1 

24.15 

1 

8 

3 

2 

1.81 

6 

2.36 

1 

0 

1.39 

1 

0 

1 

1.53 

9 

1.53 

5 

1.58 

5 

9 

3 

3 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

3.93 

7 

1.98 

3 

2.26 

7 

10 

4 

1 

(b) 

(b) 

2.29 

2 

0 

4.04 

0 

0 

1 

(b) 

(b) 

8.36 

2 

19.76 

1 

11 

4 

2 

2.60 

4 

3.96 

0 

0 

4.09 

0 

0 

0 

<b) 

(b) 

1.57 

2 

1.49 

7 

12 

4 

3 

(b) 

(b) 

1.63 

1 

0 

1.32 

7 

0 

0 

(b) 

(b) 

(b) 

(b) 

(b) 

(b) 

13 

5 

4 

2.11 

3 

1.39 

11 

0 

1.20 

7 

0 

1 

2.93 

8 

1.16 

51 

1.60 

18 

SU 

MT 


1.00 

30 

1.00 

24 

1 

1.00 

17 

1 

1 

1.00 

108 

1.00 

75 

1.00 

68 


“Number of active constraints. 
b Method failed for this problem. 

Sequence of Unconstrained Minimizations Technique. 


method gives the full set of active constraints found by 
SUMT for any of the larger problems. In each case, the best 
OC method gives the largest number of active constraints 
(method 5 for the three 60-bar truss problems and problem 
FSW.d, method 8 for IC-Wing.a, and method 13 for problem 
FSW.a). Improving the ability of these methods to drive the 
solution to contain a more appropriate number of constraints 
in the active set would probably improve the performance of 
the OC methods. 

B. Performance of Fully Utilized Design Method 

For stress constraints alone, the fully utilized design 
(method 17) produces good results for all the problems that 
have been solved. Furthermore, for simultaneous stress and 
displacement constraints, method 17 gave a satisfactory 
design for problem BarS.b (table 12), the ring (Bar60.d, 
table 24), and for problem IC-Wing.c (table 28). For prob- 
lems BarS.c (table 13) and the 10-bar truss problems 
(tables 17 to 20), the addition of a frequency constraint does 
not cause any further difficulty. However, with problem 
Bar60.f (table 24), method 17 gives a 54 percent heavier 
design than obtained with SUMT. 

To examine this last observation in more detail, consider, 
first, situations in which no stress constraints occur (or none 
become active). In this case, the design is obtained by a 


single scaling of the design variables, which normally pro- 
duces only one active constraint. It has been seen that, in the 
absence of stress constraints, the design can be significantly 
heavier. In a similar way, when more than one nonstress con- 
straint could become active (such as one displacement and 
one frequency constraint), the method 17 design is obtained 
from the fully stressed design by a single uniform rescaling 
for violated displacement or frequency constraints. This 
rescaling typically produces significantly fewer active con- 
straints than a mathematical programming technique, and 
often results in a heavier design. 

To further illustrate this aspect of FUD, recall problem 
Bar3.a (table 36), in which SUMT yields an optimal weight 
of 100 lb with six active constraints, consisting of one fre- 
quency, two displacement, and three stress constraints. The 
fully utilized design (method 17) yields an over-design 
(122 lb) with only two active constraints. Neither the fre- 
quency constraint nor any of the stress constraints becomes 
active with FUD. 

C. Performance of Hybrid Methods 

In the hybrid methods, stress constraints are treated sepa- 
rately from displacement and frequency constraints. The 
stress ratio technique is adopted for stress constraints, and 
optimality criteria techniques are followed for displacement 
and frequency constraints. Thus, when problems contain 
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stress constraints alone, the hybrid methods produce results 
identical to those of the fully utilized design. Similarly, when 
problems contain only displacement or frequency constraints, 
the hybrid methods usually produce results identical to those 
of OC methods 10, 11, and 12. We expected the hybrid meth- 
ods to work better than FUD or OC on problems in which 
both stress and nonstress constraint types appear. With these 
problems, however, the hybrid methods give unsatisfactory 
results. This may be partially explained by the generally poor 
results of methods 10, 11 and 12 in those situations. The 
hybrid methods sometimes perform somewhat better than the 
corresponding modified OC method, but sometimes worse 
(see, for example, Bar60.d and Bar60.f in table 24). Use of 
the exponential or linear form of the Lagrange multiplier 
update formulas (designated 2 and 1, respectively) may 
improve results. 

D. Convergence Characteristics 

Examination of intermediate results each time a structural 
analysis is completed during optimization of the 60-bar 
trussed ring under displacement and frequency constraints 
(Bar60.e) provides some insight into the convergence charac- 
teristics of the optimality criteria methods. Figure 13 depicts 
the weight of the design at each call to the analyzer for three 
of the OC methods, as well as for SUMT. Table 41 gives the 
total number of calls to the analyzer, as well as the total CPU 
time required to arrive at the final design when the analysis is 
run on one processor of a Cray-YMP8/8 128 running ver- 
sion 6.0 of the Unicos operating system and using version 5.0 
of the cft77 Fortran compiler. 

All four methods show oscillations initially. The SUMT 
oscillations show ever-decreasing amplitude and a slowly 
increasing mean until convergence to 392 lb is reached after 
282 analyzer calls. The OC method 16 proceeds to a point 
near 150 calls, then shows very small oscillations around 
413 lb for another 250 calls. Method 2 shows a smooth con- 
vergence to a low-weight infeasible design after 264 calls, 
then jumps to 398 lb to regain feasibility. Method 5 continues 
to oscillate throughout the optimization procedure around a 
somewhat heavier design, ending with a final weight of 
449 lb after 429 analyzer calls. The larger CPU times used by 
the OC methods were partly due to the larger number of calls 
to the analyzer, and partly due to an increased fraction of 
those calls that required constraint gradient information. 
Overall, the convergence characteristics for some of the 
optimality criteria techniques appear to be competitive with 
other nonlinear optimization methods. 

E. Number of Design Variables and Number of 
Constraints 

For an optimal design problem with n independent design 
variables, the consideration of a maximum of rt active con- 
straints is sufficient to establish the optimal point. Active 



Figure 13.— Sixty-bar trussed ring: problem Bar 60.e. 
Displacement analyzer, displacement, and frequency 
constraints. (SUMT, Sequence of Unconstrained 
Minimizations Technique; OC, optimality criteria.) 


TABLE 41 .-60-BAR TRUSSED RING: 
PROBLEM BAR60.e 
[Displacement and frequency constraints ] 


Methods 

Optimality 

criteria 

(OC) 

formulas 

Optimum 

weight, 

lb 

Number of 
structural 
analyses 

Cray-YMP 
CPU time, 
sec 

X 

X 

2 

1 

2 

397.505 

265 

136.774 

5 

2 

2 

448.536 

429 

223.209 

16 

5 

7 

412.579 

483 

251.177 

SUMT® 

391.815 

228 

85.000 


Sequence of Unconstrained Minimizations Technique. 


constraints in excess of these n , which have been termed as 
follower constraints, can be ignored without affecting the 
solution of the optimization problem (see section II.B.l). 
This feature is numerically illustrated by the three-bar truss 
problem Bar3.a (see tables 34 to 37). As mentioned earlier, of 
the 16 specified constraints, the optimal design given by 
SUMT has 6 active constraints, consisting of 1 frequency, 
2 displacement, and 3 stress constraints. The optimal weight 
is 100.042 lb, and the design variables are 1.089, 3.847, and 
3.265 in. 2 For this problem, three active constraints are suffi- 
cient to establish the optimal point. When solved by limiting 
the active constraints to three, the problem yields, via SUMT, 
a virtually identical optimal design, with a weight of 
100.074 lb and a design of 1.089, 3.847, and 3.266 in. 2 This 
indicates that for this problem of three design variables, three 
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active constraints are sufficient to fix the optimal design 
point. 

The observation that the objective function does not par- 
ticipate in the optimization process when the number of 
active constraints are equal to or greater than the number of 
design variables (see section II.B.l) is numerically illustrated 
by the three-bar truss problem Bar3.c (see tables 34, 35, 
and 38). For the problem with equal weight densities of 
0.1 Ib/in. 3 for each of the three elements, the optimal areas 
obtained are 3.299, 3.999, and 3.300 in. 2 . The optimal design 
has four active constraints. Observing the results found as the 
cost coefficients are varied over a wide range of weight den- 
sities (from 0.1 to 900 lb/in. 3 ; see table 38) shows that the 
optimal areas do not change with respect to changes in the 
cost coefficients as long as the number of active constraints 
are greater than or equal to the number of design variables 
(three). Because of the failure of SUMT in some cases for 
this problem, several other mathematical programming tech- 
niques were used in addition to SUMT and OC. It is interest- 
ing to note the performance of the optimization methods. The 
hybrid OC (method 14) and the FUD (method 17) perform 
well for all variations of the cost coefficients. As indicated, 
SUMT fails for two cases, one of which produces a negative 
weight. The IMSL routine fails for two cases, whereas the 
other SQP code fails for all six cases. The SLP code fails 
once, and the FD code fails for three cases and yields heavy 
designs for the other three cases. 

F. Efficiency of Analysis Methods 

The 60-bar trussed ring with stress and displacement con- 
straints (Bar60.d) is considered in further detail to examine 
the effects of using different analysis methods. The problem 
is solved with three OC methods and SUMT, each using the 
three analyzers discussed in section III.C: (1) the displace- 
ment method, (2) IFM, and (3) a simplified IFM. The final 
weights, number of structural analyses required, and the CPU 
times to obtain the solution on a Cray-YMP8/8128 are given 
in table 42. Differences in final weights range from less than 
1 percent for SUMT to 23 percent for method 16. The num- 
ber of structural analyses do not show significant differences, 
except for method 16. The CPU times reflect the slightly 
higher cost for using the implementation of IFM found in 
CometBoards and show a significantly lower cost for using 
the simplified IFM at each call to the analyzer. 

The significantly fewer number of analysis cycles with the 
use of the simplified integrated force method for method 16 
and the lower final weight can be explained by plotting the 
weight of each design versus the number of calls to the ana- 
lyzer (see fig. 14). Initially, the results are indistinguishable. 
Near call 40, however, both the displacement method and 
IFM cause OC to move from an infeasible design to near the 
optimal design found by SUMT. These paths then oscillated 
near the optimum until near call 200, where the solution 
moved to a heavier design. This movement may have been 


TABLE 42, -60-BAR TRUSSED RING: PROBLEM BAR60.d 


[Stress and displacement constraints.] 


Optimum weight 

Methods 


Displacement 

method 

Integrated 
force method 

Simplified integrated 
force method 

D 

B 

2 

1 

2 


431.838 

441.228 

5 

2 

2 

345.041 

360.779 

400.335 

16 

5 

7 

398.102 

398.102 

324.300 

SUMT" 

308.730 

308.729 

308.896 

Number of structural analyses required 

2 


2 

69 

69 

71 

5 


2 

805 

805 

805 

16 


7 

433 

433 

90 

SUMT* 

81 

81 

78 

Cray-YMP CPU time, sec 


1 

2 

27.342 

35.881 

15.702 


2 

2 

322.279 

421.115 

174.789 

16 

5 

7 

173.184 

224.857 

19.563 

SUMT* 

EB 


17.267 


Sequence of Unconstrained Minimizations Technique. 



Figure 14. — Sixty-bar trussed ring: problem Bar 60.d. Hybrid 
optimality criteria method 16 (Lagrange multiplier update 5 
with design variable update 7). Stress and displacement 
constraints. 
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caused by linear functional dependence among active 
constraints (see ref. 16). 

Overall, the three analysis methods appear to give similar 
solutions in most cases, but the simplified integrated force 
method is more efficient. The optimality criteria methods 
also may be more sensitive to small differences in constraint 
gradients than SUMT is. 

G. Linking of Design Variables 

The last two example problems (BarlO.c and BarlO.d) of 
the 10-bar truss shown in figure 8 (see tables 14 to 16) are 
considered to examine the effect of design variable linkage. 
Recall that each of the 10 elements is taken as a design 
variable in Problem BarlO.c, but these elements are linked to 
form only five design variables in Problem BarlO.d. The 


optimal designs and the convergence histories are given in 
tables 19 and 20, as well as in Figures 15 and 16. 

First observe that the overall optimal weight for the linked 
case (BarlO.d) is heavier than for the unlinked case by 
9.4 percent. This is expected, because the area of any given 
element within a linked design variable group is determined 
on the basis of the most critical bar element of the group. It is 
not allowed to be reduced even if that particular element’s 
contribution to the carrying of the loads is much less 
important. 

Next, consider the optimal design itself. In particular, the 
unlinked case has two element areas near their minimum 
bounds (0.103 in. 2 for element area 3, and 0, 104 in. 2 for 
element area 6). The structure obtained when these two ele- 
ments, with negligible areas, are removed becomes a deter- 
minant structure (see fig. 17). Even though the design is a 



Figure 1 5.— Ten-bar trussed ring: problem Bar 10.C. 
Displacement analyzer, displacement, and frequency 
constraints. (SUMT, Sequence of Unconstrained 
Minimizations Technique; OC, optimality criteria.) 



Figure 1 6. — T en-bar truss: problem Bar 1 0.d, Displacement 
analyzer linked design configuration displacement, and 
frequency constraints. (SUMT, Sequence of Unconstrained 
Minimizations Technique; OC, optimality criteria.) 
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Figure 17. — Eight-bar determinant structure formed by removing two elements from the ten-bar 
truss shown in figure 8. (Elements are circled, nodes are not.) 
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few pounds lighter, the potential for instability arises because 
the structure can become a mechanism when one of its 
elements fails and, therefore, may not be considered safe. 
For the linked case, the smallest element area (1.592 in. 2 ) is 
well separated from the minimum bound of 0. 1 in. 2 The opti- 
mum structure is indeterminant and does not become a 
mechanism when any one element fails. Even though the design 
with linked design variables is a few pounds heavier, it repre- 
sents a more practical and a safer design. 

In terms of computational requirements, the linked case 
(BarlO.d) requires fewer analyses than the unlinked case 
(BarlO.c) — perhaps because fewer decision (design) vari- 
ables need to be determined. In both cases, the same number 
of constraints need to be satisfied. 

The convergence characteristics of these two cases do not 
show major differences, although the linked case seems to 
produce more oscillations (see figs. 15 and 16). More numer- 
ous oscillations could be caused by an increased difficulty in 
adjusting the structure evenly because of changes to the areas 
of one element of a linked design variable group that are nec- 
essary to satisfy some constraint producing residual effects 
caused by associated changes to the areas of other elements 
of the same linked design group. 


VI. Conclusions 

When only displacement constraints are used, the OC method 
is satisfactory even for large structural systems with many 
design variables. However, the convergence behavior can 
become unpredictable if a problem has internal forces that are 
highly sensitive to design variable changes or if a fully stressed 
design (which can be considered as a heuristic OC) is mixed 
with the derivable exact stiffness OC method. 

When extended for general design application (i.e., stress 
and frequency constraints are included in addition to dis- 
placement constraints), the OC method satisfactorily pro- 
vided optimal designs for small and large structures under 
displacement and/or frequency constraints. It also is adequate 
for structures with a small number of design variables, even 
in the presence of stress constraints. However, the presence 
of large numbers of design variables and behavior constraints, 
the optimality criteria method follows a subset of constraints, 
resulting in a heavier design. 

The fully utilized design methodology was an adequate 
design tool when stress constraints dominated the design. 
Hybrid methods, as formulated, were unsatisfactory, but fur- 
ther research could be fruitful. 

The computational efficiency of the OC methods was 
similar to some mathematical programming techniques for 
displacement and frequency constraints, and the simplified 
integrated force method was more computationally efficient 
than the displacement method or IFM. All three analysis 
methods gave generally similar solutions, except, possibly, 


when the optimizer was especially sensitive to small changes 
in analysis results. 

It is both acceptable, and in some cases necessary, to limit 
the active set of constraints to no more than the number of 
design variables. Such behavior constraints need to constitute 
a linearly functionally independent set (ref. 16). 

The linking of design variables produces designs that are 
somewhat heavier, but it gives the designer added flexibility. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, May 25, 1993 


Appendix A — Symbols 


I 

areas 

Ai 

areas of individual elements 

[1M] 

diagonal matrix whose elements are 1 M,- 

a,b,c,d 

constants used in the melange form of the 
design variable update formula 

B 

force equilibrium matrix 

C 

compatibility matrix 

Q 

ith column of [C] 

Cj a 

actual value of displacement at a particular 
node or stress in a particular element 

Cja 

maximum specified permissible value of 
displacement at a particular node or 
stress in a particular element 

D 

rescaling vector 

Di 

ith component of the rescaling vector 

[D] 

matrix used in IFM for calculation of 
sensitivities 

E 

Young’s modulus 

Ei 

moduli of elasticity of individual elements 

e ij 

variables used in the illustration of a design 
variable update formula derivation 

F 

internal forces 

Ei 

ith component of the internal force vector 

IF] 

diagonal matrix used in illustrating the 
sensitivity matrix simplification 

[F] 

diagonal matrix whose elements are F, 

[F/A 2 ] 

diagonal matrix whose elements are F;M/ 2 
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ft 

fm 

fn 

fnO 

Jr 

fix) 

G 

G i 

l&l 

Sf 

Sii 

Sii 

Sis 

Sjad 

gjaf 

gjs+j 

* 

S] 

Sjix) 

8 max 

J 

[J] 

K 

X 

[M] 

[M*] 

m 

n 

njad 

njaf 


left eigenvectors of eigenvalue problem 

IFM “force mode shape” of eigenvalue 
problem 

natural frequencies of structure 
limitations on natural frequencies of structure 

right eigenvectors of eigenvalue problem 
objective function (typically weight) 
concatenated flexibility matrix 

partial derivative of the concatenated flex- 
ibility matrix with respect to the ith 
element area 

diagonal matrix whose elements are tiAAi Ej) 
constraints on frequencies 

value of /th diagonal element of flexibility 
matrix 

(iAA^Ej) (used in IFM) 

value of stress constraint corresponding to 
ith element 

constraint value for specific “active” dis- 
placement constraint 

constraint value for specific “active” fre- 
quency constraint 

jth displacement constraint 
jth active constraint 
jth constraint function 

value of largest violated displacement or 
frequency constraint 

partial derivative of the deformation coeffi- 
cient matrix with respect to the ith 
element area 

deformation coefficient matrix defined as 
([/] = m rows of [S]' T ) 

number of iterations required to reach a 
solution 

Lagrangian 

length of individual elements 
mass matrix 
IFM mass matrix 

number of inequality constraints; number of 
displacement degrees of freedom 

number of design variables; number of 
internal forces 

number of displacement constraints in 
“active” set 

number of frequency constraints in “active” 
set 


Px 

specified load in x direction 

p y 

specified load in y direction 

P 

m component load vector 

P* 

n component IFM load vector 

Po 

initial value of an update parameter, typically 0.5 

<?o 

acceleration parameters that form exponent 
of D; 

P<31 

scaling factor for ith design variable based on 
the maximum stress component associated 
with that variable 

r 

number of compatibility conditions 

Si 

partial derivative of the IFM governing 
matrix with respect to the ith element 
area 

[S] 

IFM governing matrix 

ISjg] 

diagonal matrix whose elements are (-guF^lAi 

tk 

constraint thickness at Mi iteration 

to 

initial prescribed constraint thickness 

Uj 

jth displacement component 

UjO 

displacement limitation for jth displacement 
component 

Ml 

displacement at node 1 in x direction 

«2 

displacement at node 1 in y direction 

Wj 

product of the ith bar element’s length with 
its density 

X 

m component displacement vector 

a 

step length; acceleration parameter used 
to modify update parameter, often taken 
as 1.0 

a k 

step length at the fcth iteration 

P 

acceleration parameter that is used to form 
the exponent of D, 

P 

deformation vector 

Po 

initial deformation vector 

r 

P\P 2 + 2 piP 3 + P 2 P 3 (used in the illustra- 
tion of a design variable update formula 
derivation 

SR 

effective initial deformation vector 

5 C)/(S4,) 

partial derivative with respect to the ith 
element area 

A, I 

Lagrange multipliers 

A* 

A ja 

value of Lagrange multiplier associated with 
jath “active” constraint at the fcth iteration 

X* 

vector of Lagrange multipliers associated 


with active constraints 


39 


V 

number of follower constraints; Poisson’s 
ratio 

p 

weight density 

Pi 

variables used in the illustration of a design 
variable update formulation 


design stress for ;th element 

9/o 

permissible stress for jth element 

Gu&2i> ->GLi 

stress values for each element associated 
with %i under each load condition 

X 

factor by which constraint thickness is 
reduced at each iteration 

9 

direction 

9 k 

direction vector at kth iteration 

Xi 

ith component of design variable 

Xi 

ith component of design variable at kth 
iteration 

yk+l/2 

Ai 

ith component of an intermediate design 
variable 

xt B 

ith component of the design variable lower 
bound 

yUB 
A I 

ith component of the design variable upper 
bound 

yO,k 

A'l 

ith component of the fully stressed design at 
the kt h iteration 

x) 

jth representative linked design variable 

X 

design variable 

x°* 

optimal design 

—a, opt 

fully stressed design 

r 

initial design 

(O 

circular frequency 

V/ 

gradient of the objective function 

Y/;- 

ith component of the gradient of the objec- 
tive function 

VG, 

ith column of the stress gradient (or sensi- 
tivity) matrix 

Vgj 

gradient of jth active constraint 

M 

constraint gradient (or sensitivity) matrix 

iVgja ] 

gradient of jath active constraint 

[V^Ji 

ith component of the gradient of the jath 
active constraint 

[VX] 

gradient matrix for displacement constraints 


[VXf 

approximate gradient matrix for displace- 
ment constraints 


[Ver] 

sensitivity matrix for stress constraints 


[Vflf 

approximate sensitivity matrix for stress 
constraints 


0 

zero vector 


[]-' 

inverse of a matrix 


rf 

transpose of a matrix 


ir T 

inverse transpose of a matrix 


Subscripts: 

u 

jth and jth variables 


p 

total number of stress constraints 

i 

k 

value at £th iteration 


max 

maximum 


0 

initial; permissible 


Superscripts: 



a 

approximation to sensitivity matrix 

= 

k 

value at fcth iteration (for vectors and vector 
components, e.g., <p* raised to 

the power of k (for scalars, e.g., a k , 
P\ T M ) 


L 

linked 


LB 

lower bound 


UB 

upper bound 


V 

virtual 


0 

initial 

= 

* 

constraints within active set 


Appendix B — Acronyms and Initialisms 


ANALYZ/DANLYZ Analyze/Dynamic Modes Analyze 

i 

CometBoards 

Comparative Evaluation Test Bed of 



Optimization and Analysis Routines 

I 


for the Design of Structures 

■ * 

FD 

feasible directions 

=- 

FSW 

forward-swept wing 

= 

FUD 

fully utilized design 

- 

IFM 

integrated force method 

m 

IMSL 

International Mathematical and 

- 


Statistical Library code 
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NAC 

number of active constraints 

NDV 

number of design variables 

OC 

optimality criteria method 

OC 2 

OC method using Lagrange multi- 
plier update method 1 and design 
variable update method 2 

OC 14 

OC method using Lagrange multi- 
plier update method 5 and design 
variable update method 5 

OC 17 

fully utilized design 

SLP 

sequential linear programming code 

SQP 

sequential quadratic programming 
code 

SUMT 

Sequence of Unconstrained Minimi- 
zations Technique 

VM/CMS 

Virtual Machine/Conversational 
Monitor System 


Appendix C — Derivation of Optimality 
Criteria Formulas for Displacement 
Constraints 


In this appendix, the three-bar truss shown in figure 1 is 
used as an example to derive the exponential form of the 
optimality criteria design variable update formula for dis- 
placement constraints (eq.(18)). The simplification suggested 
for the calculation of the sensitivity matrix via the integrated 
force method (see eqs. (31) and (33)) is also illustrated. 

A, Derivation of Optimality Criteria Design Update 
Formula 

First, consider the derivation of the design variable update 
formula. The original derivation of the optimality criteria for- 
mulas utilized the classical force method (refs. 3 and 4). Here 
the integrated force method of analysis is used instead to 
establish the relationship between the bar-element forces 
F and the bar-element areas X (which are taken to be the 
design variables) for any load condition (“real” or “virtual”). 
The matrices and vectors associated with the governing equa- 
tions of the integrated force method (eq.(26)) are given next 
for the three-bar truss shown in figure 1 . The force equilib- 
rium matrix is 


'Vi 

0 -£i 
2 

2 

_V 2 


2 

2 j 


The compatibility matrix and concatenated flexibility 
matrix are 




Pi 


0 


0 0 

— 0 
Pi 

0 -L 

P 3 


(C2) 


(C3) 


_ A i E i 

where Pi , and A/, E x , and ( x are the areas, Young’s 

moduli, and lengths of each of the three bar elements. 

Combining these to form the matrix S and inverting, gives 



V2 P\(p2 +P 3 ) 

pi{pi-pi) 

-V2p 3 (p 1 +p 2 ) 


-V2 Pl p 3 V2P!P 2 P3 

-P2(P1+P3) ~ 2 PiP2P3 

-V2 Pl p 3 V2 p,p 2 p 3 


(C4) 


where y = p i p 2 + 2pjp 3 + p 2 p 3 . 

For specified loads P x and P y in the x and y directions, 
respectively, and with no initial deformations, the right side 
of equation (26b) becomes 



(C5) 


Finally, dropping the zero terms causes the forces to 
become 


F =S~ l P*=- 

'fePliPl + Pi) 

-V2p,p 3 

rp r i 

Pl{P3~P\) 

_ P2(Pl + Pi) 

1 X 

p 

Y 

-V2 p 3 (p, +p 2 ) 

-V2 p,p 3 

l r yj 


Note, here, that the matrix premultiplying the load vector is 
the transpose of the matrix J seen in equations (31) and (33). 

Consider, now, the Principle of Virtual Work for the three- 
bar truss, which can be written as> 


pVpV 

r x 1 y 


u \ 

u 2 


F\F 2 F 3 v ]/j = [f/F 2 1 ’F 3 v ] GF 


(C7) 


41 


Here, u\ and u 2 are the displacements at node 1 in the x and y 
directions, respectively, the superscript v represents “virtual,” 
and the last equality comes from the discussion immediately 
following equation (26). 

For u\ y taking P x = 1 and P y = 0 and using equation (C6), 
we have 


Similar remarks can be made with respect to the displace- 
ments in the y direction. 


«2 


£12. + £22 + £32 

i4j A 2 At, 


(C14) 



_ 1 

V2 Pi ( P2 +p 3 ) 


P2(P3“Pl) 


7 

— V2p 3 (p, +p 2 )_ 


where 


(C8) 


So, substituting into equation (C7) yields 



[V2p 1 (p 2 +p 3 )p 2 (p 3 -p 1 ) -V2 p 3 (p, +P 2 )}gF 


(C 9 ) 


and 


or 


*1 V2 

e 12 = “7 PlP3^1 

*\ 7 

t 1 

e 22 = “ TT Pi (Pi + P3 ) F 1 

e 2 r 

3/2 

e 32 = —P\Plh 

£ 3 7 


d"! _ e i2 

dA i A ? 


0 = 1,2,3) 


(C15) 


(C16) 


“itf' ,3( ' ,|+ '’ 2)Fj 


For weight minimization, the objective function can be 
written as 


(CIO) 



(C17) 


Now, let 

A V2 / v 

e u = 7 Pi (P 2 + P 3 Kl 

e 21 = f-^P 2 (P3-Pl)^2 (Cll) 

f, V2 X . 

^31=“^T-— P3(Pi+P2K3 


where are the product of the bar elements’ lengths with 
their densities. 

Four cases can now be examined: when both u\ and u 2 
are positive, when they are both negative, and two cases 
when they have opposite signs. Consider the case when both 
displacements u ] and u 2 are positive. (The other cases 
become trivial extensions of this case.) When equations 
(C13), (C16), (Cl 7), and (3) are used, the stationary condi- 
tion of the Lagrangian, as given by equation (9), becomes 


Then, clearly 


Ml = 


£il + £2i + £3i 
^2 ^3 


It is now straightforward to show that 


(C12) 


“l u J0At 


■ = 0 


(/ = 1,2,3) 


(c 18) 


Solving for the areas gives 


<h*\ £ii 

dAi A? 


0 = 1 , 2 , 3 ) 


(C13) 


Aj = lyii-SL 
ijTi u J° w > 


0 = 1, 2, 3) 


(Cl 9) 
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The variables e\ } are eliminated between equations (Cl 3), 
(Cl 6), and (Cl 9) to yield the element areas, for uj bounded 
away from zero, as 


the matrix C is given by equation (C2), and the elements of 

= £. 

[ q ] are — 4 — . For these conditions, D becomes 

A f E , 


Ai = 


X 

l-i±l 


XjAJ'Vuj 
u j0 


= A 


_j = i 


(C20) 


Equation (C20) in essence represents the exponential form of 
the design variable update formula (eq. (18)), with = A> 
P= 1.0, and <?o = 2.0. Note the definition of D, in equa- 
tion (16), and of gj in equation (3). 

Having used the stationary condition of the Lagrangian 
with respect to the design variables (eq. (9)) to obtain three 
equations with which to update the three design variables, it 
is natural to look to the stationary condition of the 
Lagrangian with respect to the Lagrange multipliers 
(eq. (10)) to obtain two equations with which to update the 
two Lagrange multipliers. Although the theoretical basis for 
Lagrange multipliers updates remains somewhat unresolved, 
note that the constraint equations give 


D = — 
7 


P 1 P 2 P 3 

£ i P? 
P 1 P 2 P 3 


-V2 


(, 


Pi 


F \ P 1 P 2 P 3 

i -.2 

1 1 Pi 


-V2 


&2 P 1 P 2 P 3 ^3 PlPlPl 


y 2 
l 2 P2 

E 2 P\P2Pz 

^2 P 2 
E 2 PlPlPl 

P2 


h p\ 

^ £3 P 1 P 2 P 3 

h p\ 

£3 PIP 2 P 3 

*3 P3 2 


*p]-- 

l y 


1 -V2 1 

-3/2 2 -V 2 

1 -3/2 1 


[~ F ] 


(C23) 


where [/r] is a diagonal matrix, whose diagonal elements are 
given by p. 


L 


Pi 



= 1 


0 = 1.2) 


(C2I) 


Using J formed by transposing the matrix in equation (C6) 
and using G from equation (C3) gives 


These can then be used to provide what is essentially the 
exponential form of the Lagrange multiplier update formula 
(eq. (13)): 


A $ +1 



(C22) 


where Cj a = uj, C^ a = ujq, a = 1.0, and pq = 0.5. See 
references 3 and 4 for more details. 


B. Illustration of Simplied Sensitivity Matrix Calculation 

The simplification for the calculation of the sensitivity 
matrix suggested by equation (33) is considered next. In par- 
ticular, the product of the matrices [/] [G] [D] in equa- 
tion (31) is set to the null matrix. The justification for this 
substitution can be illustrated by obtaining the matrix D, 
which is defined immediately following equation (30), for 
the three-bar truss. The matrix S~ ] is given by equation (C4), 


[J][G][D} = 

\ V 2 pi (P2 + P 3 ) P2 (P 3 _ Pl ) -V2P3 (Pi + P3 ) 

-V 2 PjP 3 -p 2 (p|+p 3 ) -V 2 P1P3 


-L 0 

Pi 

0 — 
P 2 


0 


0 — 
P 3 


I -V2 1 

-V2 2 -V2 

I -3/2 1 


P] 


-jr [ J °o X]M-l«l 


In other words, in the calculation of the sensitivities of the 
displacement constraints, the term associated with the incre- 
ment in the forces can be set to zero. 
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Appendix D — Equations for Various 
Update Methods 

A. Lagrange Multiplier 

The equations for the Lagrange update method follow: 


jfc+l/2 k & 

Xi =*x! + 


tijad 


niad 


jad=\ 




ZjXi^O + gja^ + dxfDi 


(21a) 


X k j: l =X) a [l.O + a k Po (Cj a -C* a )} (12) ^f +, =|[^ +1/2 +^(l.O + gis j\ or X \ 


k+\!2 


(21b) 


a *+i _ 

A ,V| A ; 


-ja 


ja 


( r 

Zj£ L 

r* 

\ y J a J 


(13) 


x ^ /2 = x k (\.0 + gis ) ( 22 ) 


A*=[[ v |*] r V|*] ‘[Vi'^V/ (14) 

Here, represents the value of the Lagrange multiplier 
at the kth iteration and ja represents only those Lagrange 
multipliers associated with “active” constraints. Note that in 
equation (14) the Lagrange multiplier does not depend 
directly on its value at the previous iteration. Lagrange mul- 
tiplier update method 1 is given by equation (12), method 2 is 
given by equation (13), and methods 3, 4, and 5 are given by 
variants of equation (14). See section ILD.l for more details. 

B. Rescaling Vector and Design Variable 

Equations for the rescaling vector and design variable 
update method follow: 


ft—-* 

' Y/; 


(16) 


y * + 1 _ yk D { U P k% ) 

A I Af n 


(18) 


v ^ + l — yk 
Af Af 


10 + TF— I 0 ! -10 ) 

P do j 


(19) 


v *+i _ 

Ai 


xi 


P do 


(20) 


k 

Here, D t is the rescaling vector k and Xi is the zth compo- 
nent of the design variable at the kl h iteration. Design vari- 
able update methods 1, 2, and 3 are given by equations (18), 
(19), and (20), respectively. Method 4 uses variants of 
equations (21a) and (21b). The three hybrid methods are 
given by combining equations (18), (19), and (20), 
respectively, with equation (22). See section II. D. 2 for more 
details. 
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